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EXECUTIVE  SUMMARY 


This  report  provides  analysis  via  computer  simulation  of  pulse  stuffing 
multiplexers  and  their  interface  with  various  terminal  equipment.  The  most 
important  interface  characteri Stic  of  a  pulse  stuffing  multiplexer  is  its 
output  jitter,  characterized  by  a  continuous  change  in  the  output  frequency  as 
the  demultiplexer  attempts  to  smooth  the  effects  of  deleting  pulse  (dummy) 
bits  from  the  received  signal.  The  effect  of  this  frequency  variation  is  that 
certain  terminal  equipment  with  narrowband  clock  recovery  loops  cannot  track 
the  jitter  and  will  experience  degraded  performance  including  periodic  loss  of 
bit  synchronization. 

The  focus  of  the  analysis  is  on  the  AN/GSC-24  asynchronous  time  division 
multiplexer  (ATDM),  one  of  first  pulse  stuffing  multiplexers  used  within  the 
DCS.  Because  of  a  flexible  input/output  capability,  the  AN/GSC-24  has  found 
widespread  application  in  the  DSCS,  thus  requiring  interface  with  a  variety  of 
low,  medium  and  high  speed  terminal  and  higher  level  multiplex  equipment. 
Because  of  numerous  interface  problems  which  have  arisen  in  recent  months  with 
the  AN/GSC-24,  the  need  for  interface  testing  has  been  established  and  several 
tests  have  been  conducted  at  MilDep  test  facilities.  The  need  for  analysis 
and  computer  simulation  was  recognized  as  a  means  of  supplementing  current  and 
future  testing  and  providing  insight  into  the  performance  of  the  AN/GSC-24. 

The  simulation,  written  in  the  GASP  simulation  language,  allows  modeling 
of  a  pulse  stuffing  multiplexer,  including  structuring  of  the  TDM  frame, 
selection  of  values  for  pulse  stuffing  parameters,  and  selection  of  a  bit 
synchronization  type  as  would  be  found  in  an  interfacing  terminal  equipment. 
Specifically,  for  the  AN/GSC-24,  one  can  insert  stuff  and/or  fill  bits  into 
the  multiplexer  at  various  rates  and  remove  these  stuff /fill  bits  from  the 
demultiplexer  output.  The  simulation  output  consists  of  time  domain  waveforms 
of  the  demultiplexer's  smoothing  loop  response  or  terminating  equipment's  bit 
synchronizer  response,  from  which  jitter  performance  data  can  be  determined. 
For  those  parameters  of  interest  to  typical  AN/GSC-24  applications,  simulation 
runs  have  been  made  and  results  tabulated  or  plotted.  Performance  data  are 
given  for  low  speed  applications  of  the  AN/GSC-24  as  a  function  of  user  rate 
offsets,  multiplex  configuration  (with  and  without  fill  bits),  smoothing 
buffer  strap  settings,  bit  synchronizer  type,  and  for  tandem  connections. 
Comparisons  of  these  results  with  analytical  predictions  and  actual  test 
results  are  also  given. 

Such  simulation  results  can  greatly  reduce  the  test  and  evaluation  time 
for  AN/GSC-24  equipment  interfaces  and  for  other  pulse  stuffing  multiplexers 
(e.g.,  DRAMA,  ISTDM,  and  Adaptive  Multiplexer).  Additionally,  insight  is 
provided  into  the  reasons  for  the  observed  AN/GSC-24  performance.  One 
significant  result  of  this  new  Insight  will  be  the  improvement  of  AN/GSC-24 
multiplex  configurations,  since  jitter  performance  can  now  be  characterized 
for  any  choice  of  channel /port  rates  and  an  optimum  choice  of  multiplexer 
configurations  can  be  made  accordingly. 
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I.  INTRODUCTION 


Initial  Implementations  of  digital  transmission  In  the  DCS,  such  as  the 
Frankfurt-Koenigstuhl-Vaihlngen  (FKV )  Project  and  the  Digital  European 
Backbone  (DEB)  Stage  I,  are  using  pulse  stuffing  as  the  basic  synchronization 
technique.  The  existing  AUTOSEVOCOM  network,  which  is  largely  based  on 
50  kb/s  KY-3  operation,  is  accommodated  via  a  pulse  stuffing  interface  with 
first  level  multiplexers.  The  Digital  Communications  Subsystem  (DCSS)  of  the 
DSCS  employs  a  pulse  stuffing  multiplexer,  the  AN/6SC-24  asynchronous  TDM,  for 
all  levels  of  digital  multiplexing.  Provision  of  this  pulse  stuffing 
interface  makes  it  possible  to  use  the  medium  grade  clocks  typically  found  in 
terminal  and  transmission  equipment  and  thereby  eliminates  the  need  for  a 
separate,  highly  accurate  frequency  source.  Since  major  portions  of  the 
terrestrial  and  satellite  DCS  subsystems  will  continue  to  utilize  pulse 
stuffing  for  synchronization  in  the  foreseeable  future,  it  is  important  to 
understand  the  properties  of  pulse  stuffing  interfaces  and  anticipate  problems 
with  a  pulse  stuffing  network. 

The  pulse  stuffing  process  as  used  in  digital  multiplexing  inserts,  stuff 
bits  periodically  to  adjust  all  incoming  channels  to  a  common  bit  rate.  The 
location  of  each  stuff  bit  is  communicated  to  the  far-end  companion 
demultiplexer.  Here  the  stuff  bits  are  identified  and  deleted  to  return  each 
channel  bit  stream  to  its  original  rate.  In  the  process,  pulse  stuffing 
jitter  is  introduced,  which  results  from  the  inability  of  the  demultiplexer  to 
provide  perfectly  smoothed  clock  and  data  outputs.  This  phenomenon  can  cause 
timing  problems  in  downstream  equipments  which  must  extract  timing  from  the 
jittered  clock  and  data  signals.  A  narrowband  filter,  as  found  with  most  bit 
synchronizers,  may  not  be  able  to  track  the  high  frequency  jitter  inherent  in 
most  destuffed  clock  signals.  Such  problems  led  to  a  DCA  RaD  contract  for  a 
study  of  pulse  stuffing  interfaces  [1],  From  that  study  came  reconmendati on s 
for  specifying  a  standard  interface  between  pulse  stuffing  multiplexers  and 
bit  synchronizers. 

This  technical  note  analyzes  pulse  stuffing  multiplexers  with  the 
objective  of  giving  insight  into  the  phenomenon  of  slewing  jitter,  predicting 
performance  of  various  interfaces,  and  guiding  future  test  efforts.  The  tools 
utilized  are  theoretical  analysis  where  possible  but  computer  simulation  for 
most  applications.  A  digital  computer  simulation  of  pulse  stuffing 
multiplexers  was  previously  developed  under  the  aforementioned  RaD  contract 
[1].  This  program  has  been  expanded  by  the  authors  to  include  the  capability 
of  modeling  the  AN/GSC-24  ATOM. 

Because  of  numerous  interface  problems  experienced  with  AN/GSC-24 
operation,  the  analysis  herein  has  been  focused  on  the  AN/GSC-24.  The  main 
problem  observed  in  various  interface  tests  has  been  excessive  jitter  in 
AN/GSC-24  demultiplexer  outputs  and  resultant  degraded  performance  with 
interfacing  equipment.  The  simulation  used  herein  to  study  this  interface 
problem  allows  the  structuring  of  the  AN/GSC-24  frame  Including  stuff  and  fill 
bits,  selection  of  values  for  the  demultiplexer  smoothing  loop,  and  selection 
of  bit  synchronization  parameter  values  for  the  interfacing  equipment. 


To  run  the  AN/6SC-24  simulation,  the  stuff /fill  bit  locations  are  determined 
from  a  multiplex  plan,  which  together  with  smoothing  loop  parameters  are  input 
to  the  simulation.  The  computer  program,  which  uses  the  GASP  simulation 
language,  generates  time  domain  plots  of  the  smoothing  loop  output. 
Additionally,  the  program  allows  interface  of  the  smoothing  loop  output  with  a 
bit  synchronizer  and  can  provide  time  domain  plots  of  the  bit  synchronizer 
response.  This  latter  feature  thus  simulates  a  test  of  the  AN/GSC-24  with  any 
given  terminal  equipment.  Finally,  performance  in  a  tandemed  multiplexer 
connection  is  simulated,  thus  accounting  for  networking  applications  of  the 
AN/GSC-24. 

The  following  sections  first  describe  theoretically  the  design,  operation, 
and  performance  of  pulse  stuffing  multiplexers.  A  more  thorough  description 
of  the  AN/GSC-24  is  then  given,  including  multiplex,  demultiplex,  and 
smoothing  loop  functions.  The  simulation  program  is  generally  described, 
which,  along  with  several  appendices,  allows  the  interested  reader  to  further 
use  the  program  in  studying  the  subject  problem.  Finally,  results  of 
AN/GSC-24  simulation  runs  are  presented  for  a  number  of  configurations  of 
interest  to  DCS/DSCS  applications. 

The  choice  of  simulating  and  studying  the  AN/GSC-24  multiplexer  stemmed 
from  interface  problems  which  arose  when  the  AN/GSC-24  was  tested  in  the 
laboratory,  or  in  some  cases,  when  actually  implemented  in  the  field.  The 
advantage  of  the  simulation  described  herein  is  that  interfaces  with  pulse 
stuffing  multiplexers  planned  for  use  in  the  DCS  can  be  verified  before  final 
system  design.  Any  interface  problems  uncovered  by  the  simulation  would  allow 
changes  to  equipment  or  system  design  before  final  test  and  implementation. 

The  savings  in  time  and  money  are  potentially  great.  There  are  a  number  of 
pulse  stuffing  interfaces  that  are  planned  for  the  DCS,  in  addition  to  the 
existing  systems  mentioned  earlier.  DCS  networks  such  as  DEB,  Data 
Transmission  Network  (DTN),  and  the  ECCM  Network  all  will  utilize  pulse 
stuffing  interfaces.  The  pulse  stuffing  simulation  can  be  used  to  analyze 
these  networks  for  anticipated  interfaces,  for  both  individual  equipments  and 
for  end-to-end  connections.  Indeed,  such  an  analysis  is  planned  for  the  pulse 
stuffing  design  found  in  the  DRAMA  AN/FCC-99  multiplexer. 


II.  THEORY  OF  PULSE  STUFFING  MULTIPLEXERS 


1.  DESCRIPTION 

A  time  division  multiplexer  (TDM)  may  be  categorized  as  providing  either 
synchronous  or  non-synchronous  (or  asynchronous*)  interface  with  each  input 
data  channel.  A  synchronous  interface  specifies  that  each  source  which 
interfaces  the  TDM  must  be  timed  by  the  same  clock  which  provides  timing 
within  the  TOM.  For  this  case,  there  is  no  ambiguity  between  the  arrival  time 
of  each  data  channel  bit  and  the  time  to  multiplex  each  data  channel  bit. 
However,  if  each  data  source  is  timed  by  its  own  internal  clock,  the  interface 
is  said  to  be  non-synchronous.  For  this  case,  slight  differences  between  the 
data  source  clock  and  TDM  clock  will  arise  due  to  the  clock  inaccuracies  and 
instabilities;  the  two  clocks  will  then  wander  apart  and  eventually  will  cause 
a  timing  bit  to  be  added  or  deleted  in  the  TDM,  resulting  in  loss  of 
synchronization.  To  compensate  for  this  lack  of  synchronization  between  each 
source  and  the  TOM,  some  technique  is  required  to  convert  each  non-synchronous 
input  to  a  rate  which  is  synchronous  with  the  TDM  clock  frequency.  The 
general  approach  to  a  synchronous/ synchronous  conversion  is  illustrated  in 
Figure  1.  Each  data  source  is  timed  at  some  nominal  rate  f-j  which  will 
range  to  f-j  _^Af-j,  where  the  size  of  the  window  is  dependent  on 

clock  performance  parameters  but  is  always  much  smaller  than  f^.  The  TDM 
must  be  capable  of  accepting  each  nominal  input  rate  over  its  full  range  and 
must  provide  a  conversion  of  each  input  data  channel  to  a  rate  which  is 
synchronous  with  the  internal  clock  of  the  TDM.  Figure  1  shows  a  conversion 
from  the  nominal  data  channel  frequency  f-j  to  a  frequency  f-j '  for 
synchronous  multiplexing  in  the  TDM.  Two  commonly  applied  techniques  for 
accomplishing  this  asynchronous  to  synchronous  conversion  are  described  below. 

2.  POSITIVE  PULSE  STUFFING. 

In  the  positive  pulse  (or  bit)  stuffing  technique,  each  input  data  channel 
rate  fj  is  synchronized  to  a  TDM  channel  rate  which  i s  at  a  higher  rate  than 
the  maximum  deviation  of  the  input  data  channel  (f,  •*  Af -j ) .  Figures  2  and 
3  show  typical  block  diagrams  of  the  transmit  and  receive  sections, 
respectively,  of  a  positive  pulse  stuffing  TDM.  Input  channel  data  is  written 
into  an  elastic  buffer  at  a  rate  of  fj  +afj,  and  read  out  of  the  buffer 
at  a  rate  fj'.  Since  f-j'  is  at  a  rate  slightly  higher  than  f-j  +  Afj, 
there  is  a  tendency  to  deplete  the  buffer  contents.  Hence  the  buffer  fill  is 
monitored  and  compared  to  a  preset  threshold.  When  this  threshold  is  reached, 
a  stuff  (dummy)  bit  is  requested.  At  the  next  available  stuff  opportunity, 
the  read  clock  is  inhibited  for  a  single  clock  pulse  allowing  a  stuff  bit  to 
be  inserted  into  the  outgoing  data.  Simultaneously,  the  precise  location  of 
the  stuffed  time  slot  is  coded  into  the  overhead 


*  The  term  "asynchronous"  here  will  also  be  used  to  mean  non-synchronous. 
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ASYNCHRONOUS  TDM 


’e  2.  Positive  Bit  Stufflna  Synchronization  -  Transmit  Section 


Figure  3.  Positive  Bit  Stuffing  Synchronization.  Receive  Section 


channel  and  transmitted  to  the  receive  end  (demultiplexer)  of  the  TOM.  At  the 
receive  end,  shown  In  Figure  3,  the  received  channel  data  Is  written  Into  an 
elastic  buffer,  but  the  stuff  bits  are  inhibited  from  the  overhead  channel 
which  drives  a  clock  Inhibit  circuit.  Data  is  read  out  of  the  buffer  by  a 
smooth  clock.  Clock  smoothing  Is  achieved  by  monitoring  buffer  fill  to  drive 
a  smoothing  phase-locked  loop  (filter  plus  voltage  controlled  oscillator). 

3.  POSITIVE-NEGATIVE  PULSE  STUFFING 

In  a  positive-negative  pulse  stuffing  technique,  each  input  data  channel 
rate  f^  +afi  is  synchronized  to  a  TDM  channel  rate  f^  which  is 
identicaT  to  the  nominal  input  data  channel  rate.  Buffering  of  the  input  is 
provided  just  as  shown  in  Figure  2.  However,  now  the  buffer  contents  may 
deplete  or  spill,  since  the  buffer  is  read  out  at  a  frequency  fj,  whereas 
the  input  data  rate  may  have  either  negative  or  positive  offset  about  the 
nominal  frequency  f^.  If  the  data  source  rate  is  low,  the  buffer  will  tend 
to  deplete  and  for  this  case  stuff  bits  are  added  (positive  stuff);  if  the 
input  rate  is  high,  the  buffer  will  tend  to  spill  and  here  channel  bits  are 
subtracted  (negative  stuff)  from  the  input  rate  to  match  the  TOM  rate.  For 
the  positive-negative  bit  stuffing,  two  actions  may  occur  with  each  stuff 
opportunity:  (1)  positive  stuff  or  (2)  negative  stuff.  In  a 
positive-zero-negative  bit  stuffing  scheme,  three  actions  may  occur  with  each 
stuff  opportunity:  (1)  no  stuff,  (2)  positive  stuff,  or  (3)  negative  stuff. 
Stuff  locations  are  coded  into  an  overhead  channel  to  allow  a  receiving  end  to 
properly  remove  positive  stuff  bits  and  add  negative  stuff  bits. 

4.  BIT  COUNT  INTEGRITY  PERFORMANCE 

The  receive  end  of  the  TDM  must  properly  locate  each  stuff  bit  in  order  to 
delete  positive  stuff  bits,  add  negative  stuff  bits,  and  derive  each  original 
asynchronous  source.  Incorrect  decoding  of  a  stuff  bit  due  to  transmission 
errors  will  cause  loss  of  bit  count  integrity  in  the  channel  associated  with 
the  incorrectly  decoded  command,  since  the  demultiplexer  will  have  incorrectly 
added  or  deleted  a  bit  in  the  derived  asynchronous  channel.  In  order  to 
protect  against  potential  loss  of  synchronization,  each  stuff  location  is 
redundantly  signaled  with  a  code,  where  two  actions  (stuff /no  stuff)  are 
signalled  for  a  positive  bit  stuffing  TDM  and  where  three  actions  (no 
stuff /positive  stuff /negative  stuff)  are  signalled  for  a  positive-negative  bit 
stuffing  TDM.  After  transmission,  majority  logic  applied  to  each  code  word 
determines  the  proper  action.  The  performance  measure  is  then  given  by  the 
probability  of  incorrect  decoding  per  stuff  action,  and,  from  this 
probability,  the  time  to  loss  of  synchronization  due  to  incorrect  decoding  of 
a  stuff  may  be  obtained.  If  we  assume  independent  probability  of  error  from 
code  bit  to  code  bit,  then  the  probability  of  incorrect  decoding  for  a 
particular  stuff  action  is  given  by  the  binominal  probability  distribution. 


P(SWE)  .  Prob[ stuff  word  error]  .  N  (J)  p^l-p)*-1  (1) 

1-N  4  1 

~r 

wh ere 

N  -  code  length 

p  .  probability  of  bit  error  over  transnission  channel. 

For  small  probability  of  bit  error,  this  expression  can  be  approximated  by 

P[SWE]  -  (J)  px(l-p)N“x  (2) 

where  x  ■  (N4)/2.  Given  the  statistics  {probability  density  function)  of  the 
probability  error,  p,  one  can  find  the  mean  p,  and  then  use  p  to  find  a  mean 
P(SWE).2  Likewise,  the  variance  of  p  can  be  found  from  statistics  of  p,  and 
a  variance  of  P(SWE)  can  then  be  stated.  Knowing  the  rate  of  stuffing 
opportunities  (R,  given  in  stuff /unit  time),  the  mean  and  variance  of  the  time 
to  loss  of  BCI  due  to  stuff  word  error  can  be  calculated  as  the 


Mean  Time  to  Loss  of  BCI  -  _ 

R-PT5HF) 


Variance  of  Time  to  Loss  of  BCI  -  . 

__2 -  (4) 

R’°P(SWE) 

Application  of  these  equations  to  specific  examples  is  shown  in  reference  [2], 
5.  JITTER  PERFORMANCE. 

The  removal  of  the  stuffed  bits  at  the  demultiplexer  causes  timing  jitter 
at  the  individual  channel  output,  which  has  been  called  pulse  stuffing 
jitter.  The  actual  stuffing  rate  occurs  at  average  frequency  of  f0  stuffs 
per  second  and  is  equal  to  the  difference  between  the  multiplexer  synchronous 
rate  and  the  channel  input  rate.  In  addition,  since  there  is  a  delay  or 
waiting  time  between  the  initiation  of  a  stuff  request  and  the  time  at  which 
the  stuff  is  actually  accomplished,  the  spacing  between  stuffed  bits  varies 
about  the  average  value  of  l/f0  seconds,  resulting  in  waiting  time  jitter  at 
the  demultiplexer  output.  The  phase-locked  loop  generating  the  read  clock 
acts  to  reduce  this  jitter.  The  system  designer  must  choose  f0  large  enough 
and  the  loop  bandwidth  narrow  enough  so  that  the  jitter  is  reduced  to  an 
acceptable  value  (i.e.,  at  a  low  enough  frequency  so  that  a  bit  synchronizer 
in  downstream  equipment  is  not  affected).  There  is  also  jitter  introduced  by 
the  insertion  and  removal  of  overhead  bits  (stuff  bits,  control  information, 
frame  bits)  that  must  be  attenuated  by  the  phase-locked  loop. 


The  asynchronous  input  of  a  positive  stuffing  multiplexer  consists  of 
a  data  stream  whose  frequency  is  less  than  the  synchronous  output  data 
stream.  The  difference  between  these  two  rates  is  justified  by  discrete  one 
bit  jumps  at  the  synchronous  output.  The  timing  relationship  between  the 
asynchronous  input  and  synchronous  output  is  shown  in  Figure  4.  In  this 
figure  each  output  bit  has  a  period  which  is  4/5  the  period  of  the  input  bit 
so  that  the  instantaneous  output  rate  is  1.25  times  the  input  rate.  Every 
fifth  bit  of  the  synchronous  output  is  converted  to  a  stuff  (dummy)  bit. 
Because  of  the  rate  difference,  the  phase  of  the  output  stream  increases 
linearly  with  respect  to  the  input  stream.  When  this  phase  difference 
accumulates  to  one  bit,  the  phase  Is  reset  by  inserting  a  dummy  bit  in  place 
of  an  input  data  bit.  Pulse  stuffing  jitter  is  therefore  a  sawtooth  waveform 
with  a  peak-to-peak  value  of  one  bit  Interval  and  a  period  equal  to  the 
stuffing  rate.  The  slope  of  the  sawtooth  equals  the  frequency  offset  between 
the  input  asynchronous  rate  and  the  synchronous  output  rate. 

a.  Waiting  Time  Jitter.  For  the  jitter  waveform  to  be  the  periodic 
sawtooth  described  above,  the  stuff  bit  needs  to  be  Inserted  at  the  specific 
instant  that  the  phase  error  crosses  a  threshold.  However,  only  discrete  time 
slots  are  available  for  stuff  bits,  at  the  stuff  opportunity  rate,  fs,  and 
if  phase  threshold  crossings  do  not  occur  at  these  times  an  additional  jitter 
component  is  generated.  This  so-called  waiting  time  jitter  [3]  also  forms  a 
sawtooth  waveform  as  indicated  in  Figure  5.  The  resulting  frequency  will 
generally  be  so  low  that  the  smoothing  filter  will  have  little  effect  on 
reducing  the  peak  of  this  component  of  jitter.  The  peak  waiting  time  jitter 
may  be  found  from  Figure  6.  If  the  phase  error  crosses  the  threshold 
immediately  after  an  available  stuff  opportunity,  then  the  error  will 
accumulate  until  the  next  possible  stuff  opportunity.  The  maximum  amount  of 
phase  error  which  can  accumulate  during  this  interval  equals  the  frequency 
offset  times  the  interval  between  stuff  opportunities. 


de  -  f0  ts  bits  (5) 

where  ts  «  l/fs  sec. 

The  stuffing  ratio,  p,  is  defined  as  the  ratio  of  the  actual  stuffing  rate 
to  the  st uf f 1 ng  oppo rt un i ty  rate,  or 


The  unfiltered  peak-to-peak  waiting  time  jitter  can  then  be  expressed  as 
#*  *  ^o*s  *  p^s^s  ■  ®  (7) 
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ASYNCHRONOUS  INPUT  DATA  STREAM 


Asynchronous/Synchronous  Timing  Relationships 


AVAILABLE  PULSE  STUFFING  POSITIONS 


CORRESPOND  TO  AVAILABLE  PULSE  STUFFING  SLOTS’ 
THE  DASHED  ENVELOPE  INDICATES  THE  LOW 
FREQUENCY  COMPONENT  OF  THE  PHASE  ERROR 


Peak  Value  of  Waiting  Tine  Jitter 


Figure  7  plots  the  peak-to-peak  amplitude  of  the  filtered  waiting  time  jitter 
[3j.  Although  the  jitter  should  be  directly  proportioned  to  p,  the  smoothing 
filter  has  some  effect  on  jitter  amplitude,  particularly  at  lower  stuffing 
ratios  where  the  smoothing  filter  does  not  significantly  attenuate  the 
jitter.  At  p  -  1,  peak-to-peak  waiting  time  jitter  is  1  bit  and  its  frequency 
is  so  low  that  the  smoothing  loop  doesn't  attenuate  it.  At  p  «  1/2,  the 
peak-to-peak  jitter  equals  1/2  bit  (6  dB  below  value  at  p  ■  1);  likewise,  at  p 
»  1/3,  1/4,  and  1/5,  the  peak-to-peak  jitter  equals  1/3  (9.6  dB),  1/4  (12  dB), 
and  1/5  (14  dB)  bits.  In  general,  for  p  >  1/n,  the  peak-to-peak  jitter  equals 
1/n.  The  conclusion  reached  from  Figure  7  is  that  the  waiting  time  jitter 
peaks  for  rational  values  of  n  (with  p  -  1/n).  Hence  waiting  time  jitter  can 
be  minimized  by  selecting  a  non-rational  n  for  the  stuffing  ratio.  However, 
the  alternative  of  lowering  the  stuffing  ratio  to  minimize  jitter  would  cost 
additional  overhead  and  further  jeopardize  bit  count  integrity  performance. 


Figure  7.  Halting  Time  Jitter  vs  Stuffing  Ratio 


b.  Effect  of  Smoothing  Filter  on  Pulse  Stuffing  Jitter.  Although 
waiting  time  jitter  will  generally  occur  at  a  low  frequency  so  that  It  cannot 
be  removed  by  smoothing,  the  stuffing  jitter  component  can  be  attenuated  by 
narrowband  smoothing  loops,  as  shown  in  Figure  8. 


Figure  8.  Model  of  Smoothing  Loop  Effect  on  Jitter 

The  smoothing  loop  output  waveform  may  be  computed  by  convolving  the  loop 
impulse  response  with  the  sawtooth  phase  jitter  waveform  which  represents  the 
unsmoothed  jitter.  For  a  first  order  loop  the  closed  loop  transfer  function  is 

H(«)  -  (8) 

«c  <j« 

The  corresponding  impulse  response  is 

h(t)  -  2»fce_2,,fct  ,  (9) 


where  fc  is  the  3  dB  cutoff  frequency  of  the  smoothing  filter.  The  periodic 
sawtooth  waveform  can  be  written  as 


m(t)  m  t-n  -0.5  «  n  <  t  <  0.5  *  n  n  «  0,  J  1,  1  2,  ....  (10) 
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where  the  period  has  been  normalized  by  the  reciprocal  of  the  stuff ing^ rate, 
f0,  and  the  sawtooth  peak-to-peak  amplitude  has  been  normalized  by  360*  (1 
bit).  The  result  of  the  convolution  of  the  first  order  smoothing  loop  with 
the  periodic  sawtooth  is 

y(t)  -  /  «(«)  2»fce  da  .  (11) 

<30 

Upon  carrying  out  the  integration,  the  output  jitter  over  the  period 
-0,5  <  t  <  0.5  is 

y(t).t-  1  -e-2*V  I  ^ 

^  (  l-e-2*fc 

The  phase  jitter  y(t)  given  by  Equation  (12)  is  plotted  in  Figure  9.  Note 
that  the  smoothing  loop  with  narrower  bandwidth  provides  greater  jitter 
attenuation.  Another  important  pulse  stuffing  multiplexer  characteristic  is 
the  slew  rate,  i.e.,  rate  of  change  of  phase,  which  is  of  course  the 
definition  of  frequency,  given  by 

f (t)  =  y(t).  (13) 

Carrying  out  the  differentiation, 

f(t)  .  l-2.f_e“2"fct  I  ‘  'fc  \  (14) 

c  li^rv 


A  plot  of  f(t)  is  given  in  Figure  10  for  selected  values  of  smoothing  loop 
bandwidth.  Again  note  that  greater  slew  rate  attenutation  is  provided  by 
narrower  smoothing  loop  bandwidths. 

c.  Effect  of  Jitter  on  Bit  Synchronizers.  A  significant  problem  in 
integrating  a  pulse  stuffing  multiplexer  into  a  system  which  contains  bit 
synchronizers  is  the  effect  that  jitter  and  slew  rate  have  on  timing  recovery 
circuits.  Excessive  jitter  can  cause  a  phase-locked  loop  (PLL)  to  lose  lock 
or  slip  bits.  Since  timing  derived  by  the  PLL  is  used  to  make  bit  decisions, 
it  is  not  enough  that  the  loop  merely  maintain  lock.  It  must  also  track 
instantaneous  bit  timing  or  bit  decisions  will  be  made  with  a  timing  error. 

For  example,  with  matched  filter  demodulation  of  NRZ  data,  a  timing  error  of 
one-quarter  bit  corresponds  to  a  6  dB  loss  in  signal-to-noi se-ratio,  and  a 
timing  error  of  0.1  bit  corresponds  to  a  1.9  dB  loss  in  SNR.  The  ability  of 
the  bit  synchronizer  to  track  input  jitter  is  improved  by  increasing  the  loop 
bandwidth  but  at  the  expense  of  added  noise  bandwidth,  so  that  bit  sync 
bandwidths  are  generally  no  greater  than  necessary  for  proper  timing  recovery. 


arrutt 


Figure  9.  Stuffing  Jitter  Waveform  at  Output 
of  First  Order  Smoothing  Filter 


Figure  10.  Instantaneous  Bit  Rate  Waveform  at  the 
Output  of  a  First  Order  Smoothing  Filter 
Attributed  to  Stuffing  Jitter 


16 


In  order  to  analyze  the  effect  of  pulse  stuffing  jitter  on  bit 
synchronizers,  a  model  of  a  first  order  smoothing  loop  cascaded  with  a  first 
order  bit  synchronizer  has  been  developed,  as  shown  in  Figure  11.  The 
transfer  function  of  this  cascade  from  input  signal  to  bit  sync  phase  detecter 
output  is 

2wfc  j  2»f  (jii>) 

HU)  .  C  .  ljw)  =  C  .  (15) 

~L2wfc^JL2.k'fc-(OT 

By  partial  fractions  this  expression  can  be  expanded  to 


The  expression  is  the  linear  sum  of  the  first  order  loop  expressions.  We  can 
convolve  the  sawtooth  waveform  representing  the  jitter  input  with  this  H(w)  as 
done  before  to  arrive  at 


This  is  the  expression  for  phase  error  over  one  stuffing  jitter  period, 

-0.5  <  t  <  0.5.  This  waveform  has  a  maximum  value  at  the  extremes  (t  =  *  0.5) 
and  a  minimum  point  in-between.  The  minimum  point  is  found  by  differential’ ig 
ee(t),  setting  equal  to  zero,  and  solving  for  t.  The  minimum  point  is 
located  at 


t  =  -1/2  - 


The  peak-to-peak  phase  error  is  then 


l-e-2l,fc 


Vp  rr 


1  1  Tl-  1-e_2Wfc  <  A"k\.  fk  .  1-€"2'fc  \  ^ 

'  i_e-2wfc  i_€-2wkf„  I" J*  \  ‘  ~  1  | 


3  *  2nkfc 
k  =  3/a 


Smoothing  Loop 


Bit  Synchronizer 


FIGURE  11.  Block  Diagram  of  First  Order  Smoothing  Loop 
Cascaded  with  First  Order  Bit  Synchronizer 


Using  Equation  (19),  curves  are  shown  in  Figure  12  for  the  peak-to-peak  timing 
error  when  a  first  order  smoothing  loop  is  interfaced  with  a  first  order  bit 
synchronizer.  This  figure  presents  the  information  for  a  constant  ratio  of 
bit  synchronizer  loop  bandwidth  to  smoothing  filter  loop  bandwidth.  The 
curves  only  consider  the  effects  of  stuffing  jitter  since  the  analytical  model 
did  not  consider  waiting  time  jitter  effects.  To  illustrate  the  significance 
of  these  curves,  first  note  that  when  the  smoothing  filter  loop  has  a 
bandwidth  equal  to  the  stuffing  rate  (fc  -  f0),  the  timing  error  is 
significantly  decreased  only  when  the  bit  synchronizer  has  a  large  bandwidth 
of  approximately  10  times  the  stuffing  rate.  When  the  bit  synchronizer 
bandwidth  and  the  smoothing  filter  bandwidth  are  both  equal  to  the  stuffing 
rate  (k-1,  fc=fo)»  the  timing  error  is  approximately  30  percent  of  a  bit, 
which  corresponds  to  an  8.0  dB  loss  in  signal-to-noi se  ratio. 


PEAK  TIMING  ERROR  NORMAL tZEO  TO  PEAK  UNSMOOTHED  JITTER 


SMOOTHING  FILTER  CUTOFF  FREQUENCY  NORMALIZED 
TO  FUNDAMENTAL  FREQUENCY  OF  STUFFING  RATE 


Figure  12.  Effect  of  Pulse  Stuffing  Jitter  on  Data  Recovery 
Timing  for  a  First  Order  Bit  synchronizer 
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III.  DESCRIPTION  OF  AN/GSC-24  ASYNCHRONOUS  TDM 

Here  the  operation  of  the  AN/GSC-24  multiplexer,  demultiplexer,  and 
smoothing  loop  will  be  described.  The  AN/GSC-24  employs  positive- 
zero-negative  pulse  stuffing.  Each  channel  input  has  individual  elastic 
buffers  which  have  data  written  into  them  by  the  source  clock  and  data  read 
out  from  the  buffer  by  the  multiplexer  output  clock.  Since  the  input  and 
multiplexer  rate  elastic  buffers  are  not  synchronous,  eventually  the  buffer 
would  lose  or  gain  a  bit.  At  this  time,  a  pseudo  bit  must  be  injected  into  or 
a  data  bit  removed  from  the  data  stream  to  realign  the  phase  relation  of  the 
data  stream.  This  operation  of  inserting  or  removing  bits  must  be  controlled 
so  the  demultiplexer  knows  when  to  do  the  Inverse  operations.  This  knowledge 
is  passed  onto  the  demultiplexer  via  the  overhead  bits  which  convey  the 
stuffing  bit/data  bit  location. 

1.  MULTIPLEXER/ DEMULTIPLEXER  DESIGN 

The  GSC-24  design  is  structured  around  the  use  of  ports  for  combining 
applied  channel  inputs  into  a  single  interleaved  serial  output  data  stream.  A 
port  is  an  interval  of  time  during  which  data  from  a  particular  channel  are 
allowed  into  the  output  stream.  This  rate  is  named  the  port  rate,  and  once 
established  for  a  given  multiplexer  configuration,  remains  constant  for  all 
ports  in  use.  To  enable  the  multiplexer  to  simultaneously  accept  input 
channels  of  different  (mixed)  rates,  multiple  ports  are  assigned  to  a  single 
input  channel.  This  arrangement  is  termed  port  strapping.  For  the  GSC-24, 
the  total  number  of  assigned  ports  to  a  particular  channel  may  vary  between  1 
and  25.  The  total  number  of  assigned  ports  to  all  input  channels  of  a  given 
configuration  is  the  sum  of  all  individual  channel  values.  One  additional 
port  is  used  within  the  multiplexer  for  the  insertion  of  overhead  data  into 
the  output  bit  stream.  Thus  the  total  number  of  ports,  including  overhead, 
used  for  any  given  configuration  is  N  *1.  The  multiplexer's  output  rate  is 
the  product  of  the  total  number  of  ports  used  (N  *1)  and  the  port  rate  (Rp) . 

In  general,  a  digital  input  channel  rate  will  be  nominally  equal  to  the 
product  of  port  rate  and  the  number  of  ports  assigned  to  that  channel.  For 
cases  where  this  is  not  true,  the  multiplexer  is  capable  of  performing  a  rate 
conversion  process  termed  COARSE  RATE  CONVERSION  (CRC).  The  most  important 
benefit  resulting  from  the  multiplexer's  coarse  rate  conversion  capability  is 
that  it  enables  simultaneous  processing  of  digital  channel  input  rates  from 
mixed  rate  families  into  a  single  data  stream. 

When  the  channel  sampling  rate  for  a  specific  channel  cannot  be  configured 
to  be  within  a  tolerance  of  J-250  ppm  of  the  incoming  data  rate,  the  channel 
sampling  rate  for  that  channel  is  set  higher  than  the  incoming  data  rate.  The 
coarse  rate  conversion  circuits  then  delete  a  predetermined  number  of  gated 
clock  signals  to  effectively  lower  the  channel  sampling  rate  to  within  a 
tolerance  of  ■*  250  ppm  of  the  incoming  data  rate.  With  the  channel  sampling 
rate  then  witFin  the  J  250  ppm  tolerance,  the  stuff  command  adds  or  deletes 
gated  clocks  until  the  output  data  rate  is  synchronous  with  the  input  data 
rate.  Hence,  the  coarse  rate  conversion  circuits  are  independent  of  the 
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positive  and  negative  stuff  functions  performed  in  the  multiplexer.  The 
stuffing  of  gated  clocks  does  not  change  the  basic  bit  rate  K  •  Rp  of  the 
generated  gated  clocks  where  K  is  the  number  of  used  ports  strapped  together 
to  service  an  active  channel.  The  addition  or  deletion  of  these  gated  clocks 
is  accomplished  as  part  of  the  multiplexer's  overhead  service  capability. 

Data  leaving  the  multiplexer  is  placed  into  a  predetermined  format  as 
shown  in  Figure  13.  This  message  format  contains  data  bits  sampled  from  the 
multiplexer's  channel  data  inputs  interleaved  with  overhead  data  bits 
generated  within  the  multiplexer.  Overhead  data  are  inserted  for  two  distinct 
purposes:  (1)  to  provide  a  framing  or  synchronization  pattern  used  by  the 
receiving  demultiplexer,  and  (2)  to  inform  the  demultiplexer  of  stuffing 
actions  taken  by  the  multiplexer  as  part  of  its  input  rate  buffering 
capability. 

The  largest  grouping  in  the  output  message  format  is  the  major  frame, 
which  contains  31  complete  overhead  messages  in  addition  to  sampled  channel 
data.  Each  overhead  message  consists  of  stuffing  information  and  a  framing 
pattern  corresponding  to  one  of  the  multiplexer's  31  ports.  An  overhead  word 
Is  composed  of  one  complete  overhead  message  in  each  minor  frame.  A  minor 
frame  contains  29  data  words  and  one  complete  29-bit  overhead  word.  This  is 
the  second  largest  organized  grouping  in  the  output  message  format.  One  bit 
of  the  overhead  word  as  shown  in  Figure  14  is  located  at  the  first  bit 
position  of  each  of  the  29  data  words  in  the  minor  frame. 

Each  data  word  contains  a  sampled  channel  data  bit  for  each  port  of  the 
multiplexer  with  the  overhead  data  bit  occupying  the  first  data  word  bit 
position.  The  total  number  of  data  word  bits  is  variable,  a  function  of  the 
number  of  multiplexer  ports  in  use.  Each  multiplexer  configuration  must  use  a 
minimum  of  15  ports  giving  a  minimum  data  word  length  of  16  bits.  This  word 
length  can  range  up  to  32  bits  when  all  32  ports  are  in  use. 

When  CRC  is  utilized,  the  CRC  circuits  increment  a  word  counter  to 
generate  a  10-bit  binary  code  that  sequentially  advances  one  count  during  each 
word  until  the  end-of-scan  signal  in  word  24  of  minor  frame  count  31  occurs. 

At  this  time,  the  yord  counter  is  reset  to  zero.  The  end  of  scan  signal 
together  with  word  24  clocks  a  data  bit  into  a  5  bit  shift  register  so  an 
inhibit  signal  is  clocked  from  the  register  during  word  29.  This  signal 
inhibits  the  word  counter  and  gated  clock  control  so  a  gated  clock  cannot  be 
deleted  during  word  29  which  is  used  for  the  normal  overhead  service 
function.  Each  incremented  word  count  from  the  word  counter  is  applied  to 
Input  A  of  a  comparator.  Before  the  circuits  are  activated,  the  strapping 
switches  on  the  multiplexer  are  set  to  a  predetermined  value  related  to  the 
number  of  gated  clocks  to  be  deleted  during  the  major  frame.  This  10-bit 
binary  code  from  the  strapping  switches  is  applied  to  Input  B  of  the 
comparator  which  compares  the  10-bit  binary  word  from  the  word  counter  with 
the  10-bit  strapping  switch  code  word.  Each  time  the  binary  code  applied  to 
Input  A  is  greater  than  Input  B,  the  gated  clock  inhibit  signal  is  generated. 
In  straight  binary  logic,  the  gated  clocks  to  be  deleted  when  A  is  greater 
than  B  would  be  grouped  together.  This  would  cause  an  undesirable  series  of 
consecutively  gated  clocks.  To  prevent  this  occurrence,  the  word  counter 
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NOTES; 

1.  EACH  DATA  WORD  CONTAINS  ONE  OVERHEAO  (0/H)  BIT  PLUS  ONE  DATA 
BIT  FOR  EACH  USED  PORT.  TOTAL  NUMBER  OF  WORD  BITS  FOR  A 
GIVEN  CONFIGURATION  IS  BETWEEN  16  and  32. 

2.  EACH  MINOR  FRAME  CONTAINS  29  DATA  WORDS.  AN  0/H  BIT  FROM 
EACH  DATA  WORD  MAKES  UP  ONE  29-BIT  0/H  WORD  PER  MINOR  FRAME. 

3.  EACH  MAJOR  FRAME  CONTAINS  31  MINOR  FRAMES. 


Figure  13.  AN/GSC-24  Output  Message  Format 
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outputs  are  reversed  so  the  most  significant  bit  (MSB)  of  the  counter  Is 
applied  to  the  least  significant  bit  (LSB)  of  the  A  input  of  a  1-bit 
comparator.  This  results  in  a  homogeneous  spread  of  the  gated  clock  deletions 
over  the  major  frame. 

In  the  demultiplexer,  the  data  bits  are  normally  written  into  an  elastic 
buffer  except  for  the  stuffing  bits.  When  the  overhead  bits  identify  stuffing 
bits,  the  clock  is  inhibited,  causing  the  stuff  bit  not  to  be  written  into  the 
buffer.  If  excess  data  information  were  carried  in  these  overhead  bits,  they 
would  be  written  into  the  elastic  buffer  at  the  correct  location  of  the  bit 
stream.  This  abrupt  discontinuous  clocking  of  the  bit  stream  causes  pulse 
stuffing  jitter.  The  smoothing  loop  generates  the  clock  to  read  data  from  the 
elastic  buffer.  This  loop  contains  a  phase  lock  loop  which  attempts  to  smooth 
the  abrupt  di scontinuities  and  other  jitter.  It  is  intended  to  clock  the  data 
stream  out  of  the  buffer  at  a  stable  rate. 

2.  SMOOTHING  LOOP  DESIGN 

A  schematic  of  the  GSC-24  smoothing  loop  is  given  in  Figure  15(a).  This 
second  order  analog  phase-locked  loop  (APLL)  consists  of  a  phase  detector  that 
produces  an  error  pulse  train  whose  frequency  is  proportional  to  the  phase 
error  between  the  input  signal  and  the  voltage-controlled  multivibrator  (VCM) 
output.  The  loop  filter  has  a  low-pass  characteristic  which  averages  the 
1  phase  error  to  control  the  VCM  such  that  the  phase  error  becomes  zero.  The 

j  GSC-24  loop  filter  is  realized  with  an  active  network,  as  shown  in  Figure 

\  15(b).  Two  important  loop  characteristics  are  the  damping  factor,  n,  and  the 

i  natural  frequency  of  the  loop,  wn,  which  are  related  to  loop  parameters  by 


Kp  =  gain  of  phase  detector 
Ky  =  VCM  gain 

VCM  =  Voltage- Controlled  Multivibrator 


Figure  15.  6SC-24  Smoothing  Loop  and  Loop  Filter  Characteristics 


Another  Important  loop  parameter  Is  the  3-dB  bandwidth  [4],  which  is  given  as 
a  function  of  the  second  order  loop's  damping  factor,  n,  and  natural 
frequency,  •*: 

*3  dfl  "  “n  *  j(2n2*l)2*l  j  H  (23) 

In  applying  the  smoothing  loop  to  a  particular  user  frequency  (bit  rate), 
appropriate  choices  of  loop  parameters  are  made  by  strapping  options  in  the 
GSC-24.  Within  the  smoothing  buffer  (SB)  card  of  the  GSC-24,  the  natural 
frequency  of  the  loop  filter  is  controlled  by  switch  S5  which  allows  the 
selection  of  the  parameter  N,  while  the  damping  factor  of  the  loop  filter  is 
controlled  by  switch  S12,  which  allows  selection  of  R2  and  C2.  Resulting 
loop  parameter  values  for  selected  low  speed  data  rates  and  recommended  switch 
settings  are  given  in  Table  I.  As  part  of  the  AN/6SC-24  testing  conducted  in 
August  1980  at  the  ITF,  Ft.  Monmouth,  NJ,  certain  experimental  switch  settings 
were  tried  in  an  attempt  to  improve  interface  performance  [5].  Table  II  shows 
loop  parametric  values  which  resulted  from  optimum  switch  settings  found 
during  this  experimentation.*  Note  that  with  the  S12  switch  set  to  B',  the  3 
dB  bandwidths  of  this  smoothing  loop  for  each  data  rate  were  significantly 
reduced  as  compared  to  loop  bandwidths  for  normal  S12  settings.  As  expected, 
for  narrower  loop  bandwidths,  jitter  performance  showed  marked  improvement 
(see  Tables  IV  thru  VI). 

TABLE  I.  SMOOTHING  LOOP  PARAMETRIC  VALUES  FOR  AN/GSC-24  WITH 
RECOMMENDED  SETTINGS  (S12  SWITCH  =  C)  ON  SMOOTHING 
BUFFER  CARD 


Bit  Rate 
(b/s) 

Kp 

Ky 

ftl 

(M  ohms) 

<»F> 

N 

(k  ohms) 

W 

n 

(Hz) 

n 

“3dB 

9600 

.111 

27650 

3.3 

22 

32 

270 

1.1494 

3.414 

9.886 

4800 

.111 

27650 

3.3 

22 

64 

270 

.8127 

2.414 

4.092 

2400 

.111 

27650 

3.3 

22 

128 

270 

.5747 

1.707 

2.129 

1200 

.111 

27650 

3.3 

22 

256 

270 

.4064 

1.207 

1.146 

♦Discussions  with  Martin  Marietta  engineers  indicated  that  choice  of  non¬ 
standard  switch  settings  may  result  in  GSC-24  instability,  so  that  this 
approach  for  improving  performance  was  discarded. 
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TABLE  II.  SMOOTHING  LOOP  PARAMETRIC  VALUES  FOR  AN/GSC-24 
WITH  EXPERIMENTAL  SETTINGS  (S12  SWITCH  -  B') 

ON  SMOOTHING  BUFFER  CARD 


Bit  Rate 
(b/s) 

KP 

S 

"i 

(M  ohms) 

(»F> 

N 

R2 

(k  ohms) 

wn 

(HZ) 

n  w3dB 

(rads/sec) 

9500 

.111 

27650 

7.3 

22 

32 

T572 

1.1494 

0.2301 

1.8522 

4800 

.111 

27650 

3.3 

22 

64 

18.2 

.8127 

0.1627 

1.2863 

2400 

.111 

27650 

3.3 

22 

128 

18.2 

.5747 

0.1151 

0.9013 

1200 

.111 

27650 

3.3 

22 

256 

18.2 

.4064 

0.0814 

0.6344 

3.  INTERFACE  TEST  RESULTS 

A  brief  description  of  past  GSC-24  interface  testing  is  useful  in 
understanding  the  problem  of  slewing  jitter  and  its  effect  on  other 
equipment.  The  phase  I  (RaD)  model  of  GSC-24  employed  a  smoothing  buffer 
based  on  a  digital  first  order  loop.  A  measured  slew  rate  of  2*  in  768  bit 
times  resulted  [6],  but  since  several  sink  equipments  could  not  accept  this 
slew  rate,  modifications  were  made  to  extend  the  smoothing  interval.  However, 
this  modification  to  the  digital  smoothing  loop  also  restricted  the  allowed 
data  source  rate  offset.  For  Phase  II  (production)  units,  the  GSC-24  was 
reconfigured  and  partially  redesigned.  The  reconfiguration  consisted 
primarily  of  a  reduction  in  the  number  of  channels  (user  inputs)  from  31  to 
15,  and  the  redesign  consisted  primarily  of  a  change  from  a  digital  first 
order  loop  to  an  analog  second  order  loop  for  demultiplexer  output  smoothing. 
Figure  16  shows  the  slew  rates  determined  by  calculation  from  the  second  order 
analog  loop  equations.  In  practice,  it  is  desirable  to  measure  the  actual 
slew  rate  produced  by  a  given  smoothing  buffer.  One  successful  technique  uses 
a  delayed  trigger  oscilloscope  monitoring  the  output  timing.  One  bit,  several 
bit  times  after  the  trigger  point,  is  displayed  and  the  peak  bit-to-bit  jitter 
is  measured.  This  jitter  represents  the  phase  change  that  has  occurred  since 
the  trigger  point.  If  no  stuffing  action  has  occurred,  the  bit  edge  remains 
motionless  except  for  very  small  movement  from  random  bit-to-bit  jitter  and 
noise.  The  maximum  amount  of  jitter  is  seen  when  the  sampling  interval  on  the 
scope  occurs  during  the  steepest  slope  section  of  the  smoothing  action  (see 
Figure  9).  A  worst  case  approximation  of  the  slew  rate  can  be  found  by 


HSSiJIsEss: 


ZSSZSi 


BB5SSBB55»5b5iBBBB«B55BSS55Sp«S8SI 

liiiiaiiniuiwiuSMiiiiiiiiiiiinmniimiiiiaiHl 
■iiiiiiiiinuiiniinkNuiiiiimiiimiiiiimiiitun 


mmmm, 5SC55gS5BBBB3SjgSi 


SSESSSSSSSMBHK^n^HSSSnSSMS^MlSbiBSSSJnSSHHHH 

i 

BBMBl 

»■■■■■■!■■■■■■■— W«W»»»I>MM1WWWW»W»MM 

inainiBimmauiii^niumauiuBiiauflmiHH 

iiiuinniiiiiiuiuiimiiiimniuimiuiiiiiniip 

MMiimmnnmHHiimimmiHmwmH 


S3SS33SSSS; 


naaaniaimaij 
laaaa—aa— aiarJi 
mwaanmaiuMBi 


■niniHjamnmnninujnBa 

laaaaaaaaaamaaaKaaaaaaaa— a  *53 
■■■Bimimnnuiarinnma| 

aaaiaaapa— aaa— — R— — — — — 1 

■■■■wfijg— ■ ■ipmwBBWumwwwipiwl 

■■BBIflBBli»BIIHUB8nilBU»HBBlBBHinBBBBBB|flMDBl 

Hiuwmum— Mwimnnmmumwliimi 


i 


i 


extrapolation  of  the  steepest  slope  section  of  the  smoothing  action.  As  an 
example,  if  the  phase  change  (peak-to-peak  jitter)  were  measured  as  5%  over  a 
200  bit  interval,  then  a  100a  (2w)  change  would  occur  in  4000  bits,  giving  a 
slew  rate  of  2*  in  4000  bit  times. 

a.  HN-74  Interface  Test.  The  April-Way  1979  compatibility  test  of 
various  GFE  equipment,  including  the  MD-921/G  PSK  Modem,  CV-3034  A/D 
Converter,  AN/USC-26  Group  Data  Modem,  MD-823/G  Channel  Modem,  and  HN-74 
Interface  Device,  has  indicated  various  degrees  of  difficulty  in  interfacing 
the  GSC-24  [6].  As  an  example,  interface  tests  with  the  HN-74  operating  at 
its  standard  rate  of  1.544  Mb/s  indicated  that  the  HN-74  could  not  track  the 
frequency  variations  output  by  the  GSC-24  [7].  At  1.544  mb/s  the  GSC-24  slew 
rate  is  between  6  and  10  K  bit  times  for  a  2*  phase  change  (see  Figure  16). 

At  low  offsets  where  stuffs  occur  less  frequently  than  every  64  bits  (below  j- 
250  Hz)  the  loop  tracks  out  each  stuff  with  a  frequency  offset,  then  returns 
back  to  1.544  Mb/s.  The  average  offset  during  the  smoothing  interval  is 
roughly  150-225  Hz  with  peak  offsets  to  450  Hz.  Since  the  HN-74  can  only 
accept  offsets  of  150-250  Hz,  the  bit  synchronizer  in  the  HN-74  occasionally 
loses  phase  lock,  resulting  in  loss  of  bit  count  integrity.  This  problem  has 
been  partially  solved  by  the  use  of  coarse  rate  conversion  (CRC)  within  the 
GSC-24.  One  particular  test  configuration  utilized  100  fill  bits  per  major 
frame,  giving  a  channel  rate  of  1.566  Mb/s  in  the  multiplexer.  With  this 
increased  channel  rate,  additional  stuffing  actions  occur  during  the  nominal 
6-10  kilobit  smoothing  interval.  As  a  result  the  loop  no  longer  tracks 
individual  stuffs  but  now  goes  to  the  average  channel  rate  and  makes  small 
variations  around  it.  On  the  average,  a  stuff  action  occurs  every  70  bit 
times,  allowing  a  very  smooth  output.  Configuring  the  GSC-24  channel  in  the 
coarse  rate  mode  allowed  operation  with  the  HN-74  over  the  range  *  70  Hz 
without  degradation. 

b.  VF  Modem  Interface  Tests.  More  recently,  interface  compatibility 
tests  of  the  GSC-24  with'  several  synchronous  voice  channel  modems  have 
indicated  unsatisfactory  performance.  In  May  1979,  DCEC  issued  a  requirement 
for  GSC-24  and  voice  channel  modem  integration  testing  [5],  which  was 
subsequently  conducted  at  the  integrated  test  facility  (ITF),  SATCOMA, 

Ft.  Monmouth,  NJ.'  The  effects  of  pulse  stuffing  and  course  rate  conversion 
were  to  be  investigated.  The  multiplex  conf igurations  used  were  plans  F  and  G 
from  the  DCEC  Test  plan  [5],  as  shown  here  in  Table  III.  In  all  cases,  a 
RS-232/MIL-STD-188  interface  was  used  between  the  modems  and  GSC-24.  The  test 
configuration  used  for  each  modem  is  shown  in  Figure  17.  Modems  tested  were 
the  (1)  CODEX  9600/ V. 29,  (2)  RACAL  MILGO  MPS  9601,  and  (3)  WECO  9600. 

For  the  initial  set  of  tests,  the  modem  clock  was  not  synchronized  to  the 
GSC-24  reference  clocks,  since  this  would  be  a  normal  operational 
configuration.  Also,  normal  strap  settings  on  the  SB  card  were  used,  namely, 
strap  C  for  switch  S12  and  strap  C  for  switch  S5.  As  shown  in  Figure  17,  a 
bit  error  rate  test  set,  the  TS  3642,  was  used  to  determine  the  BER  for  each 
modem/GSC-24  configuration.  Configuration  details  and  BER  performance  are 
given  in  Tables  IV-VI  for  the  three  modems.  In  all  cases,  unsatisfactory 
performance  was  observed  for  the  three  modems  using  multiplex  plan  F  (without 
fill  bits).  The  cause  of  unsatisfactory  performance  was  a  loss  of  modem 
synchronization  due  to  GSC-24  pulse  stuffing  jitter.  This  loss  of  modem 
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synchronization  was  evident  from  bursts  of  errors,  and  correlation  of  loss  of 
synchronization  with  pulse  stuffing  events  was  evident  through  waveform 
observation  on  an  oscilloscope.  In  an  attempt  to  improve  this  observed 
performance,  other  non-standard  strap  settings  on  the  SB  board  Here  made  and 
BER  tests  conducted  with  results  as  shown  in  Tables  IV-VI.  The  strap  setting 
on  switch  S12  was  changed  from  C,  which  is  recommended  for  rates  of  1.2  to  9.6 
kb/s,  to  B  and  B',  which  are  recommended  for  higher  rates,  in  order  to 
lengthen  the  smoothing  interval  from  1000  bit  times  to  2500  bit  times.  In  so 
doing,  performance  at  7200  and  9600  b£  was  improved,  with,  in  fact,  a  zero  BER 
for  the  B‘  settings  for  all  rate  offsets  tested  and  for  all  three  modems. 

Also  for  the  rate  4800  b/s  with  switch  ST 2  again  set  with  a  B1  strap, 
performance  likewise  improved  with  zero  BER  readings,  although  it  was 
necessary  to  strap  the  S5  switch  to  the  6-12  kb/s  position  in  order  to  obtain 
a  zero  error  rate  for  the  CODEX  LSI  9600/ V. 29  modem. 

For  multiplex  plan  6,  which  used  fill  bits,  the  error  rate  did  not  appear 
to  increase  over  the  values  shown  in  Tables  IV-VI,  although  errors  still 
occurred  coincident  with  pulse  stuffing.  When  a  buffer  was  interposed  between 
the  demultiplex  channel  output  and  the  modem  data  input,  as  shown  in  Figure 
18,  error  free  operation  was  observed  for  an  interval  consistent  with  the 
clock  rate  difference  between  the  test  set  (reading  into  the  buffer)  and  the 
modem  (reading  out  of  the  buffer)  clocks.  The  buffer  compensated  for  the 
phase  slewing  action  associated  with  both  stuff  and  fill  bit  deletion  in  the 
demultiplexer.  However,  in  using  the  internal  clocks  of  the  test  set  and 
modem,  the  interval  between  buffer  resets  proved  insufficient  to  provide 
satisfactory  performance.  A  cesium  beam  standard  with  two  frequency 
synthesizers  was  then  introduced  to  provide  highly  accurate  and  stable  clock 
frequencies  to  both  the  test  set  and  modem.  As  expected,  buffer  reset  periods 
were  extended  to  operationally  useful  intervals  for  all  modems  under  test. 

Conclusions  reached  on  the  basis  of  these  observations  are: 

(1)  Pulse  stuffing  jitter  on  the  GSC-24  demultiplexer  clock  signal  causes 
frequent  loss  of  synchronization  in  each  modem's  clock  recovery  circuitry. 

(2)  For  certain  non-standard  strap  settings  within  the  smoothing  buffer 
(SB)  card  of  the  GSC-24,  acceptable  performance  was  realized  with  each  modem 
for  most  bit  rates.  This  performance  resulted  from  an  extension  of  the 
smoothing  interval  sufficient  to  allow  tracking  by  modem  bit  synchronizers. 
However,  use  of  these  non-standard  strap  settings  is  not  recommended  because 
of  unknown  performance  for  all  sets  of  configuratons,  especially  for  tandem  or 
nested  operation.  Martin-Marietta,  the  GSC-24  contractor,  was  consulted 
regarding  this  fix  to  the  problem  and  advised  against  changing  the  standard 
strap  settings. 

(3)  Use  of  a  buffer  and  a  high  quality  station  clock  provided 
satisfactory  performance,  but  represents  a  high  cost  solution. 

(4)  Approaches  other  than  those  described  above  have  been  selected  for 
interface  of  VF  modems,  including  (a)  use  of  VF  analog  channels  and  (b)  use  of 
synchronous  multiplexing,  such  as  the  LSTDM. 
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Table  III.  AN/GSC-24  Multiplex  Plan 
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TABLE  V.  RACAL  MILGO  9601  TEST  WITH  AN/GSC-24 
FOR  MULTIPLEX  PLAN  F 


Bit  Rate  Strap  Strap  Block  Length  Rate  Offset Error 
Channel  (b/s)  S12 _ S5_ _ (bits) _ (BPS) _ Rate 


1 

9600 

C 

6-1 2k 

10® 

0 

6x10“® 

B* 

6-1 2k 

10® 

4.1 

zero 

B' 

6-1 2k 

10® 

-0.9 

zero 

2 

7200 

C 

6-1 2k 

106 

0 

5xl0~® 

B* 

6-1 2k 

10® 

4.4 

zero 

3 

4800 

C 

3 -6k 

106 

0 

5X10-4 

B’ 

3 -6k 

106 

0.5 

zero 

4 

2400 

C 

1 . 5— 3k 

106 

0 

5.7x10-5 

B' 

1 . 5— 3k 

10® 

•0.3 

3.0x10-3 

TABLE  VI.  WE CO  9600  TEST  WITH  AN/GSC-24 
USING  MULTIPLEX  PLAN  F 


Channel 


1 


it  Rate 

btrap 

Strap 

Block  Lengtl 

(b/s) 

S12 

S5 

(bits) 

6-1 2k 
6-1 2k 
6-1 2k 
6-1 2k 


IV.  DESCRIPTION  OF  PULSE  STUFFING  SIMULATION 


The  purpose  of  this  simulation  is  to  compute  the  steady  state  performance 
of  «  pulse  stuffing  multiplexer.  The  system  to  be  simulated  consists  of  a 
number  of  multiplexer/demultiplexer  pairs,  referred  to  as  nodes,  cascaded  in 
series  as  would  typically  be  done  in  a  network  implementation.  A  basic 
criterion  of  the  software  simulation  is  the  capability  to  take  the  jittered 
output  of  a  node  and  make  that  the  input  to  another  node  to  permit  the 
simulation  of  cascaded  nodes. 

Figure  19  is  a  pictorial  representation  of  the  system  model  showing  two 
nodes  in  cascade  followed  by  a  bit  synchronizer.  Each  node  has  four 
frequencies  of  interest.  These  are: 

a.  The  input  frequency  (f\)  which  is  a  T1  channel  tributary  source  for 
the  first  node. 

b.  The  transport  frequency  (f2)  for  synchronous  transmission  between 
nodes,  called  the  T2  rate  frequency. 

c.  The  nominal  frequency  of  the  VCO  (f 3)  used  in  the  smoothing  loop  for 
that  node. 

d.  The  output  frequency  which  is  the  input  frequency  to  the  following 
node. 

The  parametric  studies  of  the  jitter  cascade  effect  requires  that  the 
simulation  consist  of  first  or  second  order  phaselock  loops  at  each  node 
followed  by  a  second  phaselock  loop  representing  the  bit  synchronizer. 

The  operation  of  the  elastic  buffer  in  conjunction  with  overhead  bits  to 
indicate  a  stuffing  condition  is  modeled  as  a  phase  detector  which  compares 
the  input  frequency  {f]J  with  the  node  data  transport  frequency  (f?)  and 
makes  discrete  phase  corrections  to  the  f?  input.  This  corrected  f£  clock 
rate  is  an  input  to  a  phase-locked  loop  which  models  the  smoothing  that  occurs 
in  the  system  for  reading  data  out  of  the  elastic  buffer.  The  phase  error 
sampler  is  driven  by  the  T2  rate.  The  frequency  out  cf  the  first  node  is  a 
jittered  T1  rate  which  is  an  input  to  node  2.  For  the  general  case,  the  data 
stream  continues  through  n  cascaded  nodes  and  terminates  with  a  model  of  a  bit 
synchronizer.  The  bit  synchronizer  is  another  phase-locked  loop;  the 
difference  is  that  the  parameter  of  interest  for  the  bit  synchronizer  is  the 
phase  error  {the  output  of  the  phase-locked  loop's  summing  node  or  the  input 
to  the  sampler)  which  determines  total  system  error  in  detecting  the  proper 
bit  in  the  data  stream. 

Each  node  contains  a  digital  phase-locked  loop  to  simulate  the 
demultiplexer  smoothing  loop.  The  last  node  includes  simulation  of  a  bit 
synchronizer  (receiver)  as  the  final  component.  This  node  model  is  a  linear 
model  which  is  valid  for  steady-state  performance  analysis.  Since  study  of 
the  steady-state  performance  is  the  desired  simulation  output,  the  nonlinear 
performance  of  the  phase  comparator  during  the  acquisition  phase  can  be 
ignored.  The  system  is  then  subjected  to  repeated  transient  Impacts  due  to 
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pulse  stuffing  and  fill  bit  insertions.  The  natural  simulation  model 
implementation  for  a  digital  phase  lock  loop  is  a  discrete  simulation  with 
discrete  time  inputs.  However,  this  system  is  highly  oversampled  since  the 
phase-locked  loop  bandwidth  is  very  low  compared  to  the  sample  rate.  This 
condition  permits  the  simulation  model  of  the  phase  locked  loop  to  be 
implemented  as  a  continuous  system  defined  via  differential  equations  to 
provide  faster  solutions  on  the  computer  with  identical  dynamic  response.  The 
variable  step  size  Runga-Kutta  integration  algorithm  permits  relatively  large 
time  increments  for  integration  after  the  transient  from  each  pulse  stuffing 
or  fill  bit  insertion  activity. 

The  dynamic  equations  are  table-driven  to  permit  generalized  routines  for 
first  or  second  order  nodes  that  are  accessed  based  on  the  number  of  nodes 
being  simulated  for  any  one  simulation  run.  Discrete  events  include  the 
opportunity  to  advance  or  retard  the  node  phase  and  the  periodic  collection  of 
phase  and  frequency  statistics  for  analysis  and  plotting.  Minimum  and  maximum 
phase  peaks  are  detected  which  trigger  a  state-event  for  collection  of 
statistics  on  peak  phase  variation. 

To  represent  each  hardware  module  set  (node),  two  items  are  required: 

(1)  A  discrete-time  simulation  of  the  periodic  phase  correction 
opportunities  with  a  decision  criteria  to  determine  the  phase  correction  time. 

(2)  A  continuous  (either  state  equation  or  difference  equation) 
simulation  of  the  smoothing  loop's  effect  on  the  phase  relationship  of  the 
jittered  signal  passing  through  the  smoothing  loop  with  respect  to  the  stable 
input  signal. 

Since  we  are  modeling  a  linear  system  (with  additive  elements),  a 
perturbational  model  of  phase  jitter  is  the  preferred  simulation  approach  to 
guarantee  computational  accuracy.  This  is  accomplish  by  setting  the  T1  rate 
input  frequency  equal  to  zero. 

The  cascade  effect  is  modeled  by  using  the  first  node's  output  as  the 
second  node's  input  and  repeating  for  the  number  of  simulated  nodes.  Figure 
20  is  a  pictoral  representation  of  the  perturbational  model  showing  two 
cascaded  nodes. 

To  complete  the  overall  model,  it  is  necessary  to  model  the  bit 
synchronizer  as  the  final  element  of  these  cascaded  nodes.  The  bit 
synchronizer  is  a  phase-locked  loop  whose  parameter  of  interest  is  the  phase 
error,  i.e.,  the  output  of  the  phase-locked  loop's  summing  node. 

The  simulation  objectives  are  to  obtain  answers  to  the  pulse  stuffing, 
fill  bit  problems  of  various  multiplexers  as  functons  of  multiple 
configurations,  input  rates,  output  rate  offsets,  and  internal  strapping 
options  in  smoothing  loops.  This  is  accomplished  by  measuring  peak  jitter, 
rms  jitter,  and  frequency  slew  rate  versus  time  from  plots  generated  by  the 
GASP  simulation  program.  These  plots  show  various  phase  locked  loop 
parameters  plotted  as  a  function  of  time. 
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V.  RESULTS  AND  CONCLUSIONS  FROM  AN/GSC-24  SIMULATION 


The  simulated  configurations  were  extracted  from  multiplex  plans  found  in 
reference  [5].  Actual  hardware  tests  were  conducted  in  August  1979  at  the 
Integrated  Test  Facility  (ITF),  Ft.  Monmouth,  NO,  as  reported  in  Chapter  III, 
Section  3b  of  this  report.  Specifically,  multiplex  plans  F  and  G  of  reference 
[5]  were  selected  for  simulation.  A  comparison  of  simulation  results  with 
actual  results  could  then  be  accomplished. 

For  several  parametric  variations,  the  demultiplexer's  smoothing  loop 
performance  has  been  characterized  by  simulation.  This  characterization  took 
the  form  of  (1)  peak-to-peak  jitter  measured  from  the  time  response  waveforms 
output  from  the  smoothing  loop  or  (2)  slew  rate  defined  as  the  maximum  rate  of 
phase  change  and  measured  from  the  loop  output  waveforms  as  a  a0  degree  (or 
radian)  phase  change  in  a  aT  number  of  bit  times  (or  seconds).  Since  a  90* 
phase  change  would  result  in  significantly  degraded  performance,  the 
measurement  of  slew  rate  was  based  on  the  minimum  number  of  bit  times  observed 
for  a  90*  phase  change.  To  facilitate  comparison  with  the  test  results  and 
other  specifications,  the  resulting  slew  rate  numbers  have  been  extrapolated 
to  2*  radians  by  multiplying  the  number  of  bit  times  measured  times  four. 

Several  computer  runs  were  completed,  mostly  at  channel  rates  of  4800  and 
9600  b/s,  to  comply  with  the  selected  multiplex  configurations.  For  a  given 
computer  run,  depending  on  the  parameter  being  varied,  output  data 
corresponding  to  a  fraction  of  a  second  up  to  several  seconds  were  examined. 
Plots  of  smoothing  loop  phase  error  versus  time  were  used  in  determining 
peak-to-peak  jitter  and  slew  rate.  Results  of  these  measurements  for  the 
parameters  under  study  are  presented  below. 

1.  EFFECT  OF  STUFF  VERSUS  FILL  BITS 

Table  VII  indicates  results  of  testing  multiplex  plan  F  (with  no  fill 
bits)  versus  multiplex  plan  G  (with  fill  bits)  for  rates  of  4800  and 
9600  b/s.  Nominal  (recommended  in  GSC-24  Operations  Manual  [8]  parametric 
values  were  used  for  the  smoothing  loop  bandwidth  and  damping  factor.  Note 
that  multiplex  configurations  which  employed  fill  bits  (for  rate  conversion) 
were  found  to  have  better  slew  rate  characteristics.  Intuitively,  a  larger 
number  of  fill  or  stuff  bits  should  result  in  smaller  phase  changes  required 
to  keep  the  smoothing  loop  output  frequency  about  some  average  (and  desired) 
value.  A  fewer  number  of  fill  or  stuff  bits  would  result  in  larger  phase 
transients  about  the  desired  frequency.  Given  this  intuitive  reasoning,  slew 
rate  performance  should  improve  as  the  input  frequency  is  offset  by  greater 
amounts,  since  larger  offsets  will  result  in  a  higher  stuffing  ratio  for  a 
positive-negative  pulse  stuffing  multiplexer.  Examination  of  Figure  22  shows 
this  effect,  where  improved  slew  rate  performance  is  observed  for  J250  ppm 
offsets  as  compared  to  ■•10  ppm  offsets  for  up  to  four  nodes  in  cascade. 
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2.  EFFECT  OF  POSITIVE  VERSUS  NEGATIVE  STUFFING 

Table  VIII  gives  results  of  several  simulation  runs  aimed  at  showing  the 
effects  of  positive  versus  negative  pulse  stuffing  on  the  slew  rate 
characteristic.  These  data  were  collected  by  offsetting  the  input  frequency 
by  equal  amounts  but  in  opposite  directions.  For  a  positive  offset,  the 
GSC-24  will  use  only  positive  stuffing;  for  negative  offsets,  only  negative 
stuffing  will  result.  These  two  effects  are  due  to  the  positive-negative 
nature  of  GSC-24  pulse  stuffing.  Since  the  smoothing  loop  is  subject  to  the 
same  absolute  magnitude  of  phase  error  for  equal  but  opposite  offsets,  the 
slew  rate  characteristic  should  be  unaffected  by  the  sign  of  the  offsets. 
Examination  of  Table  VIII  indicates  that  indeed  the  slew  rate  is  not 
materially  affected  by  change  in  the  sign  of  the  offset. 

3.  EFFECT  OF  DAMPING  FACTOR 

Table  IX  shows  the  effect  of  changing  the  damping  factor  in  the  GSC-24 
smoothing  loop  on  the  slew  rate  characteristic.  Earlier  tests  of  the  GSC-24 
interfaced  with  VF  modems,  conducted  at  the  ITF,  Ft.  Monmouth,  NJ,  indicated 
improved  performance  with  judicious  choice  of  damping  factor.  Specifically, 
strap  S12  of  the  smoothing  loop  board  controlled  the  damping  factor.  Use  of 
the  B'  strap  resulted  in  pronounced  improvement  as  shown  in  Tables  IV  through 
VI.  Likewise,  simulation  results  shown  in  Table  IX  indicate  significant 
improvement  with  a  choice  of  the  B'  strap. 

These  data  also  illuminate  the  relationship  between  the  GSC-24  slew  rate 
and  typical  VF  modem  bit  synchronizers.  It  is  known  that  VF  modems  can  accept 
typically  up  to  J  w  radian  phase  change  in  10,000  bit  times.  A  look  at  Table 
IX  shows  that  in  every  case  a  choice  of  the  C  strap  for  damping  factor  results 
in  an  excessive  slew  rate.  Tables  IV  through  VI  reflect  this  same  conclusion, 
since  all  modems  tested  with  the  C  strap  resulted  in  a  non-zero  bit  error 
rate.  Conversely,  a  choice  of  the  B'  strap  resulted  in  acceptable  slew  rate, 
as  shown  in  Table  IX  and  as  indicated  by  the  zero  error  rates  of  Tables  IV 
through  VI. 

4.  EFFECT  OF  TANDEMING 

Figures  21  and  22  illustrate  the  effect  on  jitter  performance  of  tandeming 
GSC-24  multiplexer/demultiplexer  pairs  (nodes).  For  this  illustration,  a  9600 
b/s  data  rate  is  used  from  multiplex  plan  G  (see  Table  III),  which  was 
selected  so  that  both  fill  and  stuff  bits  would  be  utilized.  Figure  21  plots 
peak-to-peak  jitter  versus  tandemed  nodes  for  up  to  six  nodes,  while  Figure  22 
plots  slew  rate  versus  tandemed  nodes  for  up  to  eight  nodes.  In  obtaining 
these  data,  the  actual  data  rate  used  was  offset  from  9600  b/s  by  up  to  250 
ppm.  The  curves  labeled  *250  ppm  resulted  from  offsetting  each  nodal  clock 
(f 2 •  the  transport  frequency  —  see  Figure  19)  by  <250  ppm.  Likewise,  those 
curves  labeled  -*10  ppm  or  _^50  ppm  resulted  from  alternating  the  sign  of  the 
offset  with  eacF  succeeding  node.  In  this  way,  the  effects  of  rate  offsets, 
in  both  magnitude  and  sign,  can  be  observed  with  tandemed  nodes.  First  note 
that  slew  rate  and  jitter  performance  degrades  with  increased  tandeming  for 
250  ppm  offsets,  regardless  of  sign.  Results  for  10  ppm  offsets  actually  show 


slight  improvement  in  slew  rate  performance  with  increased  tandeming.  These 
results  indicate  that  jitter  and  slew  rate  performance  is  not  strongly 
dependent  on  the  number  of  tandemed  nodes.  This  same  weak  dependence  on 
tandeming  has  been  observed  with  testing  of  other  pulse  stuffing  TDM  equipment 
[9]. 

5.  EFFECT  ON  BIT  SYNCHRONIZERS 

Curves  showing  the  interaction  of  a  bit  synchronizer  with  a  jittered  input 
signal  are  shown  in  Figures  23  and  24.  These  curves  show  peak-to-peak  phase 
error  and  maximum  slew  rate  measured  at  the  output  of  the  bit  synchronizer  for 
the  GSC-24  second  order  smoothing  loop  interfaced  with  a  second  order  bit 
synchronizer.  For  these  curves  a  nominal  input  frequency  of  9600  Hz  with  an 
offset  of  0.1a  was  utilized.  The  two  figures  present  the  information  for  a 
constant  ratio  of  bit  synchronizer  loop  bandwidth  to  smoothing  filter  loop 
bandwidth,  with  the  ratio  denoted  by  k.  These  curves  show  that  jitter 
performance  is  improved  by  increasing  the  loop  bandwidth  of  the  bit 
synchronizer.  For  example,  from  Figure  23,  when  the  smoothing  loop  has  a 
bandwidth  equal  to  one-tenth  the  stuffing  rate,  the  phase  error  ranges  from 
150*  for  k=0.093  down  to  4.5*  for  k=0.93. 


TABLE  VII.  EFFECTS  OF  STUFF  VS.  FILL  BITS  ON  SLEW  RATE 


Test  Bit 

Run  Rate 

1  4800 

2  4800 

3  9600 

4.  9600 


Mux$  Fill  Bits /  Oamping  Factor  Loop  Band-  Slew 

Plan  Major  Frame  Strap  (S12) _ width  (Hz)  Rate 

F  0  C  25.7  2«  rad / 

1450  bits 

G  833  C  25.7  2w  rad / 

7600  bits 

F  0  C  62.1  2*  rad / 

3200  bits 

G  767  C  62.1  2tr  rad/ 

6250  bits 


TABLE  VIII.  EFFECT  OF  POSITIVE  VERSUS  NEGATIVE  STUFFING  ON  SLEW  RATE 


Test 

Run 

Bit 

Rate 

Mux  [5] 
Plan 

Fill  Bits/ 
Major  Frame 

Damping 

Factor 

Positive  or 
Negative  Stuff 

Slew 

Rate 

2 

4800 

G 

833 

C 

♦ 

2i  rad/ 
7600  bits 

9 

4800 

G 

833 

C 

- 

2*  rad/ 
7600  bits 

6 

4800 

G 

833 

B‘ 

+ 

2*  rad/ 
33700  bits 

10 

4800 

G 

833 

B 1 

— 

2*  rad/ 
32500  bits 

1 

4800 

F 

0 

C 

+ 

2»  rad/ 
1450  bits 

11 

4800 

F 

0 

C 

— 

2*  rad/ 
1460  bits 

5 

4800 

F 

0 

B' 

+ 

2*  rad/ 
12200  bits 

12 

4800 

F 

0 

B* 

— 

2«  rad/ 
11700  bits 

TABLE  IX.  EFFECTS  OF  DAMPING  FACTOR  ON  SLEW  RATE 


Test  Bit  Mux  [5]  Fill  Bits/  Damping  Factor  Loop  Band-  Slew 


Run 

Rate 

Plan 

Major  Frame 

Strap  ( S12 ) 

width  (Hz) 

Rate 

1 

4800 

F 

0 

C 

25.7 

2*  rad/ 
1450  bits 

5 

4800 

F 

0 

B’ 

8.1 

2w  rad/ 
12200  bits 

2 

4800 

G 

833 

C 

25.7 

2*  rad/ 
7600  bits 

6 

4800 

G 

833 

B1 

8.1 

2*  rad/ 
33700  bits 

3 

9600 

F 

0 

C 

62.1 

Zi  rad/ 
3200  bits 

7 

9600 

F 

0 

B1 

11.6 

2*  rid / 
86500  bits 

4 

9600 

G 

767 

C 

62.1 

2*  rad/ 
6250  bits 

8 

9600 

G 

767 

B* 

11.6 

2*  rad / 
53760  bits 
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Bit  rate  =  9600  b/s 
Stuff  rate  =  9.6  Hz 
Multiplex  Plan  G 
Damping  Factor  B' 
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Figure  21.  Effecting  of  Tandeminq  on  Pulse  Stuffing  Jitter  for  AN/GSC-24 


FIGURE  22.  EFFECTING  OF  TANDEMING  ON  PULSE  STUFFING  SLEW  RATE  FOR  AN/GSC-24 
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Figure  24.  Effect  on  Pulse  Stuffing  Slew  Rate  on  Clock  Recovery 
for  a  Second  Order  Bit  Synchronizer 
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APPENDIX  A 


USER'S  GUIDE 


1.  INTRODUCTION 

User  familiarity  with  the  ITEL  AS-5  Batch  and  Tine  Sharing 
Option  (TSO)  Operating  Systems  including  the  IXLBUB  Text  Editor 
is  assumed. 


2.  AVAILABILITY 


These  programs  which  include  the  Stuff/Fill  hit  and  GASP 
Simulations  are  available  for  use  on  the  DCEC  ITEL  AS-5 
computer  system. 

3.  STUFF /FILL  BIT  PROGRAM  DATA  CARD  FORMAT 


The  following  data  cards  are  required  to  use  STUFF ILL: 


Card  #1:  Channel  Information  Card 

Variable  Column  position 

NCHANS  not  specified 

NCHAN  not  specified 


Format 

free 

free 


Card  *2:  Stuff/fill  card 

Variable  Column  positions 

NFLBT  1  to  70  by  5 


Forma  t 
1415 


Card  *3  Fill  bits/channel  Card 

Variable  Column  positions 

NFL BTC  1  to  70  by  5 


Pornat 

1415 


Card  # 4 

Stuf fs/channel  Card 

Variable 

Column  positions 

NSTUFF 

1  to  70  by  5 

Card  #5 

Wait  tine  jitter  Card 

Variable 

Column  positions 

NNTJ 

1  to  70  by  5 

Format 

1415 

Format 

1415 


Card  16 


Port  and  word  Card 


A-l 


Variable 

Col ana  positions 

Fornat 

BPOBTS 

not  specified 

free 

HHOBQS 

not  specified 

free 

Card  #7 

port  to  channel  Card 

Variable 

Colusa  positions 

Fornat 

IPTCHV 

1  to  70  by  5 

1415 

4.  STUFF IL  USAGE 


Back  ran,  corresponding  to  the  generation  of  the  stuff/fill  bit 
locations  in  the  GSC-24's  coabined  bit  strean,  requires  all  seven 
inpat  cards  described  above. 

To  run  STOFFIL,  enter  HYLBUfi  and  execute  the  PDS  nenber 
B3893.LIB (HIHGO).  A  typical  session  is  listed  in  Appendix  B. 

The  input  data  for  this  prograa  typically  is  saved  in  a  partitioned 
data  set  (PDS)  whose  nenber  nane  is  specified  by  the  user. 

All  Fortran  source  prograns  are  listed  in  Appendix  D. 


APPENDIX  B 


SAMPLE  DIALOG  TO  OS  I  STOFF/FILL  PBOGBAN 

This  dialog  is  valid  only  after  the  user  has  success  folly  logged 
onto  TSO  and  has  received  the  BEAD!  response  froa  the  coaputer.  For 
this  ezaaple,  the  user  naae  is  coraack,  the  charge  nuaber  is 
11352000,  and  the  user's  badge  auaber  is  B3893.  The  stuff/fill  bit 
locations  in  the  GSC24*s  coabined  data  streaa  Mill  be  determined  as 
a  function  of  the  inputs  described  in  the  User's  Guide.  Oser 
responses  follow  the  guestion  narks. 

M 

COMMAND?  execute  froa  *r3893. lib  (nihgo) '  clear 

YOLOBE  IS  TSO  NIC  1 

3  CHAB  JOB  IDENT?  sfb 

BON  TIBE(SBC)?  60 

DATA  DSN  BEBBEB  BARE?  gsc24a 

GSC24  .  F8XX.  DATA  HORBEB?  7 

DATA  DSN  BEBBEB  NAHB? 

//B3893SPB  JOB  ( 1 1352000, 2G05, 60, 8,1 92) ,' JOHN  J  CCHHACK', 

//  NOTIFT=B3893,  HSGCLASS=Q 
//PBINT1A  EXEC  PGH=0PBINT, HEGIOR=60K 
//SYS PBINT  DD  STSOOT=A 

//SIS0T1  DD  DSN=  B3893.  DATA  (GSC24  A)  ,  DISP=S  HR 
//PBOGBON  EXEC  PGB=  GSC24  ,  HEGION=  19  2K 
//STEPLIB  DD  DSN=  B3893.  LOAD,  DISP=S  HR 
//FT04F00  1  DD  SYSO0T=A 

//FT05F001  DD  DSN=  B3893.  DATA  (GSC24  A)  ,  DISP-SHfi, LABEL-  (,  ,  ,  IN  ) 

// FT 06F00  1  DD  SYSOOT=A 

//FT08F001  DD  DSN=H3893.  GSC24.  F807.  DATA,  DISP=SHH 
// FT09F00 1  DD  DSN=S TEMPI ,DISPs(NEV, DELETE) ,SPACE=  (CYL,1)  , 

//  DCB=  ( BECF  R=  YS)  ,0NIT=SISDA 

//FT10F001  DD  DSN=STEBP2,DISP=  (NEB, DELETE)  ,SPACE=(CTL,  1)  , 

//  DCB=(HECFH=YS)  ,OHIT=SISDA 

//FT11F001  DD  DSN=6TEHP3,  DISP=  (NEB,  DELETE),  SPACE=  (CTL,1), 

//  DCB=  (RECFB=YS)  ,DNIT=SISDA 

//FT  12F00 1  DD  DSN=6TEBP4, DISP* (N EH, DELETE) ,SPACE=  (CYL,1), 

//  DCB=(BECFB=YS)  ,0NIT=SYSDA 
//PBINT1B.  EXEC  PGB=UPHINT,  HEGION=60K 
//SXSPBINT  DD  SYSOOT=  A 

//SIS0T1  DO  DS  N=  B3893.  GSC24.  F807  .  DAT  A,  DISP=SHR 

// 

BON  JOB(Y  OB  N)  ?  y 
JOB  S0BRITTED 
EXEC  END 


appendix  c 


SARPLE  DIALOG  TO  OS  E  THE  GASP  PBOGBAf! 

This  dialog  is  valid  oal y  aftsr  the  usee  has  successfully  logged 
onto  TSO  and  has  received  the  BEAD!  response  froa  the  cosputer.  for 
this  e  maple  the  user  *s  a aae  is  Coraack,  the  charge  nuaber  is 
113S2000,  and  the  user's  badge  nuaber  is  B3893.  The  first  BYLBOB 
exec  file  sets  up  the  input  data  as  required  by  the  GASP 
subroutines.  The  second  BYLB UB  exec  file  executes  the  progran  in 
the  background. 

A.  FIRST  BYLBOB  EXEC  FILE 

v 

COBH  AN  D?  execute  fron  '  r3893 .  lib  (gaspnOd)  '  clear 
VOLUHE  IS  TSOHK1 
DSN?  data 

HERB  EH  NAME?  gasp96ab 

J  J  CORHACK  24  12  27  1979  1000000000000(1 

J  J  CORHACK  24  06  10  1980  100000000000011 

2SSYNC  11  0.  520. 

LONE B  VALUE?  -360. 

UPPER  VALUE?  360. 

2SSYNC  1  1  -360.  360. 

3OO0TPUT  1  1  0.  360. 

LONER  VALUE?  -  180. 

UPPER  VALUE?  180. 

300UT  PUT  1  1  -180.  180. 

4  0  1.0E-7  1.0E-8  1.0E-8  1.0E-4  -1.0 

ABSOLUTE  LOCAL  TRUNCATION  ERROR? 

RELATIVE  ERROR? 

BIN  STEP  SIZE? 

NAX  STEP  SIZE? 

4  0  1.0E-7  1.0E-8  1.0E-8  1.0E-4  -1.0 

1  7813  5 

FRAHE  START  NUHBER?  1 
FRAME  STOP  HUH  BE  R?  500 

1  500  5 

9. 58907E- 3  1 

NBPO?  25 

9.  58907E  -3  25 

4  2 

NODE  TYPE?  2 

4  2 

9.58907E-3  25 
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.  yj;  -v 


9.6E-3 

9.6001E-3 

9.5999E-3 

9.6E-3 

9.5E-3 

T 1  BATE?  9.6e-3 
9.6E-3 
9.6001E-3 
9.5999E-3 
9.  6E-3 
9.6E-3 
9.6E-3 


09.58907B-3 

09.58907E-3 

09.58907E-3 

09.58907E-3 


25 

0  9.6E-3 
0  9.6E-3 
0  9.6E-3 
0  9.6E-3 


9.6E-3  0  9.6E-3 

T2  RATE?  9.6001e-3 
9.6001E-3  0  9.6E-3 


9.6001E-3  0  9.6E-3 

T2  RATE?  9.6002e-3 
9.6002E-3  0  9.6E-3 


9.5999E-3  0  9.6E-3 

T2  RATE?  9.6 00Ue-3 
9.6004E-3  0  9.6E-3 


9.6E-3  0  9.6E-3 

T2  RATE?  9.6 006e-3 
9.6006E-3  0  9.6E-3 


0.0  90.0 
0.0  90.0 
0.0  90.0 
0.0  90.0 

INITIAL  PHASE  OF  OUTPUT?  45.0 


0.0 

45.0 

0.0 

45.0 

0.0 

45.0 

0.0 

45.0 

1 .8293E-3 

.230701 

1.8293E-3 

.230701 

1.8293E-3 

.230701 

1.  8  293E-3 

.230701 

FN?  1 . 8294e-3 
1.8  294E-3 
1 .8294E-3 
1.8294E-3 
1.8294E-3 
ZETA?  .230701 
1.8294E-3  .230701 

1.8294E-3  .230701 

1.8294E-3  .230701 
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1.8294E-3  .  230701 

PREFIX  -  DATA  :  HEHBER  -  GASP96AB 
NONE  =  NO  SAVE,  STOP  =  PHOGBAH  END 
HEHBER  NAHE?  gzsp96ba 
HEHBER  GZSP96BA  SAVED  IN  R3893.DATA 
DSN?  stop 
EXEC  END 

B.  SECOND  HYLBOB  EXEC  FILE 

w 

COHH AND?  execute  fro*  •  r3893  .lib  (gaspgo)  •  clear 

VOLOHE  IS  TSOMK1 

3  CHAR  IDENTIFIER?  run 

RON  TIHE  (SEC)?  120 

OOTPUT  LI  NES  (K)  ?  5 

PD S  OB  SEQ?  pds 

DATA  INPOT  HEHBER  NAHE?  gzsp96ba 
GSC24.F8XX. DATA  FILE  NOHBER?  7 
LIST  INPOT  DATA  SETS?  y 

#  OF  PLOTS?  1 

#  OF  CYL  FOB  TEflP  DSN?  3 
BIT  RATE  (INTEGER)?  9600 
PLOT  EVERY  XTH  POINT?  25 
POS  OR  NEG  STOFF?  pos 
ROBE  RONS?  y 

//R3893R0N  JOB  (1 1352000, 2G05,  120,  5,  150)  , 'JOHN  J  COHN ACK* 
//  HSGCLA  SS= Q, NOTIFY=B389  3 
//PRINTAO  EXEC  PGH=0PBINT,REGION=60K 
//SYSPRI  NT  DD  SYSOOT=A 

//SYSOT1  QD  DS  N=  B3893.  DAT  A  (GZS  P96B  A)  ,  DISP=SHR 
//PBINTBO  EXEC  PGB=OPBINT,  BEGION=60K 
//SYSPBINT  DD  SYSOOT=A 

//SYS0T1  DD  DS  N=  B3893.GS  C24.  F807  .  DAT  A,  DISP=SHR 
//GO96001  EXEC  PGH=G  AS  PS  IH,  HEGION=  150K,  T  IH  E=  (  2,  0) 
//STEPLIB  DD  DSN=R3893. GASP. LOAD, DISP=SHB 
//FT03F001  DD  * 

25,-1 

/* 

// FT 0 4 FO 0 1  DD  DSN =H 3893. G SC24.F 807 .DATA ,DI SP= SHR 
//FT05F00 1  DD  DSN=R3893.  DATA  (  GZSP96BA)  ,DISP-SHR, 

//  LABEL=  (,  ,  ,  IN) 

//FT06F001  DD  SYSOOT=  A 

//FT11F001  DD  DS  N=  8G PLOT  1,  DISP=  (NER,  DELETE)  , SPACE- (CYL,  3) 
//  DCB*(RECFH=VS)  ,ONIT=SYSDA 
// 

RON  JOB(Y  OR  N)  ?  y 
JOB  SO SHITTED 
EXEC  END 
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APPENDIX  D 


STUFF/FI LL  S OURCE  LISTING 

C  DETERMINES  WHICH  DATA  WORD  CONTAINS  FILL  BITS  AND  WHERE  IN  THE 

C  DATA  WORD  THE  FILL  BITS  ARE  LOCATED 

DIMENSION  IFALAG  (3 1 )  ,  IWORD(32)  ,  MULT(  10)  ,NBFT  (10)  ,NBIT  (1  0)  ,NB  YT(  1  0) 
DIMENSION  IPTCHN  (31)  ,NFLBTC(31)  ,  M  B  FT  (  5)  ,MBIT(5)  ,  MBYT  (5)  ,NFLBX(31) 
DIMENSION  N  PTPCH  (  3 1)  ,NSTUFP(31)  ,NZEBO(14)  ,NDELPH(14)  ,NCL0CK(14) 
INTEGER  NWTJ(31),AREV,  AAREV 
DATA  MOLT/1,2,  4,  8,  16,32,64,128  ,256,512/ 

DATA  NWT  J/3  1*0/,  NS  TOFF/ 31*0/,  NPTPCR/31*0/ 

C  NCHANS  IS  THE  TOTAL  NUMBER  OF  CHANNELS 

C  NCHAN  IS  THE  CHANNEL  UNDER  INVESTIGATION 

READ  (5,*,END=999)  NCHANS, NCHAN 

C  NFL BTC  IS  THE  ARRAY  THAT  TELLS  ROW  MANY  FIL1  BITS  APE  IN  EACH 

C  CHANNEL  VIA  A  SPECIAL  CODE  WORD  FOR  EACH  CHANNEL 

READ(  5,  120)  (NFLBT(IAB)  ,IAB  =  1  ,  NCHANS) 

READ  (5,120)  (NFLBTC(IAB)  ,  IAB=1,  NCHANS) 

READ(  5,  120)  (NSTUFF(IAB)  ,IAB=1  , NCHANS) 

READ  (5  ,120)  (NWT  J  (IAB)  ,  IAB=1,  NCHANS) 

120  FORM  AT  (  141 5) 

C  N PORTS  IS  THE  TOTAL  NUMBER  OF  PORTS 

READ  (5,*)  NPOSXS  ,  NSOBDS 
C  WRITE(  13,510) 

510  FORMAT  (•  FORTRAN  FILE  13*) 

C  WR ITE ( 1 4, 6 10) 

610  FORMAT ( *  FORTRAN  FILE  14') 

WRITE  (6,190)  NPORTS,  NCHANS, NWORDS,  NSTUFF(NCHAN)  ,NWTJ  (NCHAN)  , NCHAN 
190  FORMAT  (13, '  PORTS,  ‘,12,'  CHANNELS,  ',13,'  WORDS  AMD  ',15, 

1*  CLOCKS/STUFF  AND  ',15,*  CLOCKS/WAIT  TIME  FOE  CHAN', 13) 

IF (NPORTS. LT. 15. OR. NPORTS. GT.31)  GO  TO  999 

ICLOCK=0 

LTOT  AL  =0 

ITOTAL=0 

K  FL  A  G=  0 

MMWORD=0 

NPT S=0 

READ  (5,120)  (IPTCHN  (I  AD)  ,  I  AD=  1 ,  NPORTS) 

DO  130  IAC  =  1, NCHANS 

WRITE  (6, 140)  I  AC,  NFLBTC(IAC)  ,  NFLBT  (I  AC) 

140  FORM  AT  (IX,  'CHANNEL  ',12,'  WHOSE  FILL  BIT  CODE  IS  ',14,'  HAS  ',13, 

1  •  FILL  BITS') 

130  CONTINUE 

DO  150  IAC=1, NPORTS 
WRITE(6,  160)  IAC, IPTCHN  (IAC) 

160  FORMAT  (1  X, 'PORT  ',12,'  IS  ASSOCIATED  WITH  CHANNEL  '  ,12) 

150  CONTINUE 

DO  205  I AC=1, NCHANS 
DO  205  IACC*  1, NPORTS 

IF  (IPTCHN (I ACC)  .  EQ.  IAC)  NPTPCH(IAC)  =NPTPCH(IAC>  ♦  1 
205  CONTINUE 
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DO  210  IAC= 1 , NCHANS 
MBITE(6,  215)  IAC,8PTPCH  (IAC) 

215  FORMAT  {*  CHANNEL*  ,  13 ,  •  HAS  ',12,'  POETS') 

210  CONTINUE 

LCLOCK=  8PTPCH  (NCHAN) 

C  88  IS  THE  DATA  WORD  POINTER 

HKNT=0 

DO  535  NN=1,  NWORDS 
8=88 

AN  =  FLO  AT  (8) 

DO  555  KNX= 1,10 
NBFT(KNX)  =0 
8BIT  (KNX)=Q 
555  NBTT  (KNX)  =0 

1=0 

10  I=IO 

NDI V=8/2 

NBIT (I)  =8-  2*NDIV 
N=  NDI  V 

IF(NDIV.EQ.O)  GOTO  20 
GOTO  10 

C  CALCULATES  THE  BINARY  VALUE  OF  A  IE  THE  RHOLE  NUMBER  PART 

20  DO  30  JJJJ=1,10 

30  NBYT(JJJJ) =N8IT ( 1 1- JJJJ) 

ANFRCT=  (AN-AINT  (AN)  ) 

DO  40  KKKKK=1,  10 
A=ANFRCT*2. 

IF  (A.  GE.  1.0)  GOTO  35 
NBFT  (KKKKK) =0 
AN  FRCT=A 
GOTO  4  0 

C  CALCULATES  THE  BINARY  VALUE  OF  A  IE  THE  DECIflAL  EART 

35  NBFT (KKKKK) =1 

AN  PRCT=A-  1 . 

40  CONTINUE 

ARE V=0 

C  CALULATES  THE  REVERSE  BINARY  VALUE  OF  A 

DO  67  NNDEX= 1 , 10 

67  AB  EV=AREV*8  B YT (NNDEX)  *  M  ULT  (NNDEX) 

DO  151  IAA=1,32 
IF (  IAA .  GT.  3  1)  GO  TO  145 
IFALAG  (I  AA )  =0 
145  IHORD(  IAA)  =0 

151  CONTINUE 

C  IAE  IS  THE  CHANNEL  POINTER 

IAE=NCHA8 

C  AT  EACH  DATA  IORE  CHECK  FOR  FILL  BIT  INSERTIONS 

VHITE{6,  40  1)  NN,AREV,IAB,NFLBTC  (IAE)  ,  NBIT,  NBTT 
401  FORHAT  (•  NN=  »  ,  I3,2X,  '  AR  EV=  •  ,  14,  2Xf  »  CHA8N  EL  ',12,'  FILL  BIT  CODE 
1*,  14,  '  BIN  =  ',10I1,'  REV  BIN=' ,1011) 

IF  ( ARE  V.  GT.  NFLBTC(IAE)  )  IFALAG  (IAE)=1 
VUMB1  =NPTPCH  (IAE) 


D-2 


-  > 


DO  235  LADD= 1 , NUM81 
IF  (LADD.  GT.  1 )  GO  TO  240 
IF(IFALAG(IAB)  .EQ.1)  I  WORD  (LADD)  =1 
IF  (IFALAG  (I  AE)  .  EQ.  1 )  MKNT=MKNT-»1 
IF(  IFALAG(IAE)  .EQ.O)  I  WORD  (LADD)  =0 
GO  TO  235 

240  IWORD  (LADD)  =0 

235  CONTINUE 

IFALAG  (I  AE)  =0 

C  ****************** ********************************** 

C  A  DATA  WORD  HAS  BEEN  GENERATED 

C  ************************** ************************** 

NUMB2=  NPTPCH  (NCHAN) 

DO  110  KOA=1,NUMB2 
LTOTAL=  LTOT  ALf 1 

IF (MOD (LTOTAL, N  STUFF(NCHAN))  .NE.O)  GOTO  110 
IWOBD  (KOA)=IWORD(KOA)  *1 
110  CONTINUE 

NU  MB3=  NPTPCH (NCHAN) 

WRITE  (6,200)  NN,  (IWORD  (IAF)  ,  IAF=  1,  NUHB3) 

200  FORMAT  (•  DATA  WORD  ',13,'  IS  *,32l1) 

C  ICNT  =  0  WHEN  FIRST  FILL  BIT  PER  WORD 

C  ICNT  =  1  WHEN  NOT  FIRST  FILL  BIT  PER  WORD 

ICNT=0 

C  DO  LOOP  SAMPLES  THRU  ALL  BITS  OF  DATA  WORD 

NU  HB4=  NPTPCH  (  NCH  AN) 

DO  50  JJC=1 , NUMB4 

IF (IWORD  (JJC) . LT. 1)  GOTO  180 

C  ICLOCK  IS  POINTER  TO  ARRAY  CONTAINING  FIIL  BIT  INTERVALS 

ICLOCK=ICLOCK*1 

C  IF  (  ICNT  .  EQ  .  0)  WRITE  (6, 405)  ICNT, JJC , LCLOCK,MMWORD 

4  05  FORMAT {’  ICNT=  » ,12 ,5X , » JJC= • , 14, 5X  , ' L CLOCK= • , 13, 5X, » MM  WORD  =  • , 13) 

IF  (  ICNT.  EQ.  0)  N  CLOCK  (ICLOC  K)  =J  JC  +  (NPTPCH  (NCHAN)  -LCLCCK)  ♦  MRVOBD 
IF  (ICNT. EQ.O)  RMWOBD=0 

IF  (KFLAG.EQ.O.  AND.  ICNT.  EQ.O)  10  FS  ET=NCLOCK  ( ICLOCK) 

C  IF ( ICNT. EQ . 1)  HR ITE (6,500)  ICNT , JJC , LCLOC K 

500  FORMAT  ('  ICNT=  *  ,15, 5X,  •  JJC='  ,  15,  5X,  'LCLOCK  =  '  ,  15) 

IF  (ICNT.  EQ.  1)  NCLOCK  (ICLOC  K)  =JJC-LCLOCK 
NDELPH  (ICLOCK)  =IWORD(JJC) 

IF  (ICLOCK.  NE.  1 )  GO  TO  2222 
ISTORE=N  DEL  PH  ( ICLOCK) 

C2222  WRITE (6  ,2333)  ICLOCK, NDELPH( ICLOCK)  , ICLOCK,  NCLOCK (ICLOCK) ,JJC,IWOR 
C  ID(JJC) 

2  33  3  FORMAT  (•  NDELPH  (',  12  ,')  =  ',  12,  5X ,  •  NCLOCK  (•  ,  12,'  )  ='  ,  1 2, 5X  ,  •  I  WORD  (•  ,1 
12,')=*,  12) 

2222  KF  LAG=  1 

LCLOC  K=  JJC 
ICNT= 1 

230  IF  (ICLOCK.  LT.  14)  GO  TO  55 

220  WR  ITE ( 9)  (NCLOCK  (I AH)  ,1  AH=  1  ,IC LOC  K) 

C  WRITE  (13  ,505)  (NCLOCK  { IAH)  ,  IAH=  1,  ICLOCK) 

50  5  FORM  AT  (1415) 


WBITE(IO)  (NDELPH  (ZAH)  ,IAH=  1, ICLOCK) 

C  WRITE(  14,  606)  (NDELPH  (IAH)  ,IAH  =  1  , ICLOCK) 

606  FORMAT  (1415) 

80  FORM  AT {  141 5) 

DO  250  I AZ= 1 , ICLOCK 
250  ITOTAL=ITOTAL*NCLOCK  (IAZ) 

NPTS=NPTS+ 14 
ICLOCK  =  0 
GO  TO  55 

180  IF(JJC.EQ.N  PTPCH (NCHAN)  .AND . IC NT .EQ.O )  MRWORD=  HHWORD+NPT  PCH  (NCHAN) 

55  IF  (NN.  NE.  NWORDS)  GO  TO  50 

IF (JJC .N E.NPTPCH (NCHAN ) )  GO  TO  50 
C  LAST  WORD  CASE 

225  ICLOCK= ICLOCK  ♦  1 

C  IF  (IC ‘IT.  EQ.O)  8  BIT  E  (6 , 4  10)  ICNT,  LCLOCK,  HHWOR  D,  10  FSET 

410  FORHAT  (•  ICNT  =  *  , 13 , 5X, • LC LOC  K=  •  ,13  ,5  X,  ‘HHWORD=  *  ,13 ,5X,  •  IOFS  ET=  • ,  13 

1) 

IF  (ICNT.  EQ.O)  N CLOCK  (ICLOC  K)  =(NPTPCH  (NCHAN)  -ICLOCK)  ♦  HHWORD+5*NPTPC 
IH(NCHAN)  ♦IOFS  ET 

C  IF(  ICNT.  EQ.  1)  WRITE  (6,415)  ICNT,  JJC  ,LCLOCK,IOPSET 

415  FORHAT  (•  ICNT=  •  ,  13  ,5  X  ,  »  JJC=  •  ,  13,  5X  ,  •  LCLOCK=' ,  13,  5X,  •  I0FSET=*  ,13) 

IF  (  ICNT.  EQ.  1)  NCLOCK  (ICLOCK)  =JJC-LCLOC  K*5*NPTPCH(  NCHAN)  ♦IOFSET 
NDELPH  (ICLOCK)  =I80RD(JJC) 

IF  (NDELPH(ICLOCK).  EQ.O)  N  DELPH  ( ICLOCK)  =ISTORE 
C  KRITE(6,  2432)  ICLOCK, NDELPH  (ICLOCK)  , ICLOCK,  NCICCK  (ICLOCK)  ,  JJC,  IHOR 

C  ID  (JJC) 

2432  FORM  AT  (•  LAST  HOR  D  .  .  .  .  NDELPH  ( *  ,12  ,  •  )  =  •  ,1 2  ,5  X  ,  •  NC  LOC  K  (  •  ,  12,  • )  =  •  ,  12, 
15X,  'I  H  CRD  (  •  ,  12  ,  •)  =  •  ,  12) 

WHITE  (9)  (NCLOCK  (IAH )  ,  IAH=  1  ,IC  LOCK) 

C  WFITE(  13,  505)  (NCLOCK  (IAH)  ,IAH  =  1  , ICLOCK) 

WRITE  (10)  (  NDELPH  (IAH)  ,IAH=1,  ICLOCK) 

C  WHITE{  14,606)  (NDELPH  (IAH)  ,IAH  =  1  , ICLOCK) 

DO  260  I AZ=1, ICLOCK 
260  ITOTAL  =  ITOTAL*NCLOCK  (IAZ) 

NPTS=NPTS* ICLOCK 
GO  TO  535 
50  CONTINUE' 

535  CONTINUE 

WRITE  (6  ,305)  HKNT 

305  FORH AT ( 1  TOTAL  NUMBER  OF  FILL  BITS  IS  *,I5) 

WRITE  (6, 245)  NPTS 
245  FORHAT  (•  NPTS  ='  ,15) 

END  FILE  9 
REWIND  9 
END  FILE  10 
REWIND  10 
JCOU  NT=0 
ICOU  NT=0 
NTOT=0 
NPTSS=0 
IF  LAG=0 

DO  1000  IA=  1, NPTS,  1 4 


0-4 


1045 

1025 

1035 
1  040 

1030 

1050 

C 

82 

1015 

1000 

85 

81 

83 

270 


NDIFF=NPTS-IA 

NDIFF1  =  NDI FF* 1 

IF(HDIFF.GE.  13)  BD1FF1  =  14 

READ  (9)  (NFLBT  (IAD)  ,  IAD=  1,  NDIFF1) 

READ (10)  ( IFALAG (IAD) ,IAD=1 ,NDIFP1) 

IF  (IA.  NE.  1.  AND.  NDIFF.  GE.  13)  NL  ASTP=NL  AST 
IF (N  DIFF . GE .  13)  NLA ST=NF LB T { 14 ) 

DO  101^  r.B=  1 , NDIFF1 
NTOTP-  f.  icX 

IF  (NFLBT  (IB).  NE.  1)  NTOT=0 
NTOT=N  TOT+IFALAG ( IB) 

IF  (NFLBT  (IB).  EQ.  1)  I FLAG=  1 
IF(NFLBT(IB)  .EQ.1)  GO  TO  1015 
ICOUNT=ICOUNT*  1 
NC  LOCK  (ICOD  NT)  =  NFLBT  (IB) 

IF(NDIFF.LE.  13.  AND.  IB.  EQ.  NDIFF1)  GO  TO  1045 

IF  ( ICO  UN  T.  L  T.  14)  GO  TO  1025 

WRITE  (11)  (NCLOCK  (IAD)  ,  IAD=  1,  ICOUNT) 

ICODNT=0 

IF  (IB. EQ.1  )  GO  TO  1035 

IBM  1  =  IB-  1 

NPREV=  NFLBT  (IBM1) 

GOTO  1040 

IF  (I  A.  NE.  1  )  N  PP.EV=  NL  AST  P 
IF  (  IA.  EQ.  1)  NPREV  =  1 
IF  (IFLAG.EQ.  0)  GOTO  1030 

IF(NFLBT(IB)  .NE.  1. AND. NPREV. EQ.1  )  GO  TO  1030 
GOTO  1015 

IF(  IFLAG.EQ. 0)  JCO  ON  T= JCOUNT+  1 
IF (I FLAG. EQ.O)  NDELPH (JCOU NT) = NTOT 
IF  (IFLAG. EQ. 1 )  NDELPH (JCOU NT)  =  NTOT P 
IF  (IFLAG.EQ.  1 )  JCOUNT= JCOUNT* 1 
IF  (IFLAG.EQ. t )  NDELPH (JCOUNT) = NTOT 
IFLA  G=0 

IF(NDIFF.LE. 1 3 . AN D .1 B . EQ . NDIFF  1 )  GO  TO  1050 

IF  (JCOUNT.  LT.  14)  GOTO  1015 

WRITE  (12)  (NDELPH  (IAD)  ,  I A  D=  1 ,  J  COUNT) 

NPTSS=NPTSS* JCOUNT 

WRITE  (6  ,82  )  NPTSS,  JCOUNT,  NDIFF  1,  I  A,  IB 
FORBAT('  NPTSS  =  •  ,I5,5X,  •  JCOUNT  =  ',I5, 

1 5X , '  NDIFF1  =• ,I5,5X, 'I  A  = •, 13, 5X, • IB  = •, 13) 
JCOUNT=0 
CONTINUE 
CONTINUE 

WRITE  (6,85)  NPTSS 
FORBAT  (  •  NPTSS  =•  ,15) 

WRITE  (8,81  )  NPTSS 
FCRB  AT  ( 1 10) 

IF  (NPTSS. EQ.O)  WRIT  E  (6 , 8  3) 

FQ  RB  AT  ( '  PROGRAB  TERBINA  TED  SINCE  NPTSS=0‘) 

WRITE  (6  ,270)  ITOTAL 

FORBAT  ('  THE  TOTAL  NUMBER  OF  CLOCK  PERIODS  IS*, I7) 
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IP  (NPTSS. EQ.O)  GO  TO  999 

END  FILE  11 

BEHIND  11 

IND  FILE  12 

BEHIND  12 

DO  300  IAH=1, NPTSS, 14 

NDIPF=NPTSS-IAH 

NDIFFl =NDIFF*1 

IF (NDIFF.GE. 13)  NDIFF1=14 

READ{  1  1)  (NCLOCK  (IAD)  ,IAD=1  ,NDIFF1 ) 

HRITE  (8,80  )  (NCLOCK  (IAD)  ,  IAD=  1,  NDIFF  1) 
300  CONTINUE 

DO  545  IHOL=1,14 
545  NZERO(IHOL)  =NHTJ(NCHAN) 

C  HRITE  (6,585)  NZERO(14) 

585  FORM  AT  ( '  N Z EBO  ( 1 4)  =•  ,  1 5) 

C  HRITE (6  ,80)  (  NZERO(IADDD)  ,  IADDD=  1,  1  4) 

DO  340  IAH= 1 , NPTSS  ,14 
NDIFF=NPTSS~IAH 
NDIFF1  =  NDIPF*1 
IFJNDIFF.GE.  13)  NDIFF1=14 
C  HRITE  (6,84)  NDIFF1 

84  FO  RM  AT  ( '  NDIPP1=',I5) 

HRITE  (8  ,80)  (NZERO(IAD)  ,IAD=  1,  NDIFFl ) 

340  CONTINUE 

DO  320  IAH=1, NPTSS, 1  4 
NDI FF=  NPTSS -I  AH 
NDIFF 1 =N  DIFF ♦ 1 
IF  (NDIFP.GE.13)  NDIFF1=14 
C  HRITE(6,  800)  NDIFFl 

800  FORMAT  (•  AT  READ  FILE  12.. .NDIFFl  =  •  ,  15) 
READ  (12)  (NDELPH  (IAD)  ,IAD=1 , NDIFFl ) 

HRITE (8, 80)  (NDELPH (IAD) , IAE=1  , NDIFFl ) 
320  CONTINUE 

HRITE  (4,81)  NPTS 
REHIND  9 
BEHIND  10 

DO  400  I  AH=  1,NPTS,  14 

NDIFF=  NPTS -I  AH 

NDI FP 1=N  DIFF ♦ 1 

IF  (NDIFF.GE.  13)  NDIFF1  =  14 

READ  (9)  (N  CLOCK  ( IAD)  ,IAD-=1  ,  NDIFF  1 ) 

HRITE  (4, 80)  (NCLOCK  (IAD)  ,  I  A  D=  1,  N  DIFF  1) 
400  CONTINUE 

DO  425  I AH= 1, NPTS, 14 

NDIFP=NPTS-IAH 

NDIFF1  =  NDIFF«-1 

IF (N  DIFF. GE  .  13)  NDIFF1=14 

HRITE  (4  ,80  )  (NZERO  (IAD)  ,  IAD=  1,  NDIFFl) 

425  CONTINUE 

DO  420  IAH=1,  NPTS, 14 
NDIFF=NPTS- IAH 
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HDIFP1=HDIFF*1 
IF (NDIFP.GE.  13)  NDIFF1=14 
READ  (10)  (HDELPH  (I  AD)  ,  IA£^  1,  HDIFF1) 
MBITE(4,  80)  (HDELPH  (IAD)  ,IAD=1  ,NDIFF1 ) 
420  COMTINUE 
999  STOP 
END 
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APPENDIX  E 


GASP  OS  EB  SOOBCE  LISTING 


C*****pOLSE  JITTER  STODT ************************* *********************** 

C 

PURPOSE 

SIMULATE  CASCADED  NODES  FOR  ANALYSIS  OF  WAITING  TIRE  JITTER 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


DESCRIPTION  OP  PARAMETERS 

DH 2(  )  -  DATA  TRANSFER  FREQUENCY  (BAD/SEC)  AT  EACH  NODE. 

PERTOREATION  FROM  BASIC  CHANNEL  PREQOENCY. 

DW3  {  )  -  VCO  FREQUENCY  (RAD/SEC)  AT  EACH  NODE. 

PERTURBATION  FROM  EASIC  CHANNEL  FREQUENCY. 

GAIN  (  )-  GAIN  OF  PHASELOCK  LOOP. 

TAU (  )  -  TIME  CONSTANT  FOR  SECOND  ORDER  PHASELOCK  LOOP. 
PHERR  {  )  — S  AMPLED  PHASE  ERROR  (RAC). 

TS(  )  -  TIME  BETWEEN  SAMPLES  FOR  PHASELOCK  LOOP.  PERIOD  OF 

SYNC  FREQUENCY  FOR  EACH  NODE. 

NFD  (  )  -  FIXED  DELAY  IN  PHASELOCK  LOOP  CORRECTION. 

DELAY  WITHIN  FRAME  IN  BITS. 

CHFREQ  -  CHANNEL  FREQUENCY. 

NBPO  -  NUMBER  OF  BITS  PER  OPPORTUNITY. 

NFCUB  -  CURRENT  FRAME  BEING  PROCESSED. 

NFSTRT  -  STARTING  FRAME  FOR  DATA  COLLECTION. 

NFSTP  -  STOPPING  FRAME  FOR  DATA  COLLECTION. 

DTP  -  TIME  INCREMENT  USED  TO  DETERMINE  FSLEW. 

NSpF  -  number  of  samples  per  frame. 

IDECSN  -  FLAG  FOR  DECISION  TIME.  RESET  BY  SAVE  EVENT. 

SYSTEM  EVENTS 

SAVE  -  SAVE  TIME  FOE  DATA. 

EVENT  OCCURS  EVERY  SAVE  TIME  DTSAV . 

EVENT  CODE  1. 

DCBCT  -  OPPORTUNITY  TO  DECIDE  TO  CORRECT  PHASE  OF  SYNC 
FREQUENCY  FOB  NODE.  OCCURS  EVERY  NBPF  BITS. 

EVENT  CODE  2. 

CHCTN  -  TIME  TO  CORRECT  PHASE  OF  SYNC  FREQUENCY  FOR  NODE. 
THIS  EVENT  GENERATED  FROM  DCRCT  EVENT. 

EVENT  CODE  3. 

SEVNT  -  STATE  EVENT.  DETECTION  OF  OUTPUT  PHASE  MINIMUM 
OR  MAXIMUM.  EVENT  CODE  4. 

FSAVE  -  SAVE  FREQUENCY  TO  DETERMINE  FSLEW. 

EVENT  CODE  5. 


FILES 

FILE 


1  -  EVENT  FILE. 

ATRIE(I)  -  TIME  OF  EVENT. 
ATBIB  (2)  -  EVENT  CODE. 

ATRIB  (3)  -  NODE  INVOLVED  IN 


EVENT  ACTIVITY. 


DATA  COLLECTION 


E-l 


COLCT. 
1  - 
2  - 
TINSA. 
1  - 
2  - 

3  - 

4  - 

5  - 

6  - 

7  - 

8  - 


PHHAX. 

PHNIN. 

PHOOT. 

FNODE. 

FSLEH. 

PHERRN. 

PHTRK. 

PHSHCE. 

FTRK. 

FSLEHT. 


HIS  TO  GRANS 


1  -  PHHAX. 

2  -  PHNIN. 

3  -  PHOOT. 

4  -  PHERRN. 

5  -  PHTRK. 

6  -  PHSHCE. 

7  -  FNODE . 

8  -  FS  LEV. 
PLOTS 

PLOT  1. 

3  -  PHOOT. 

4  -  PHSRCE 

5  -  PHERRN 

6  -  PHTRK. 
PLOT  2. 

1  -  F  NODE. 

2  -  F SL EH. 


PHOOT. 

PHSRCE. 

PHERRN. 

PHTRK. 


RENARKS 

THIS 

DATA 


HODEL  IS  A  PERT  OR  BAT  ION  AL  HODEL. 

IS  COLLECTED  AT  LAST  NODE  IN  CASCADE. 


SOB  PROGRANS  REQUIRED 

I  NT LC  ,E  V  NTS  ,ST  ATE, PHASE 


GASP  SUBPROGRAHS 


C  *** ****************************************************** ************** 

C 

C  NAHE  BOX  JIT 

C 

C  PROGRAH  DATE  780417 

C 

REAL*  16  QF2 
CONNON  QSET  ( 200) 

CONNON  /GCON1/  ATRIB(25)  ,  J  EV NT ,  N F A,  N  FE{  100)  ,  NLE  (100)  ,MSTOP,  NCRDR  ,N 
1NAPO,NNAPT,NNATR,NNFIL,NNQ (100)  , NNTRY,NPRNT,PPARH(5  0,  4)  ,T NOW , TTBEG 
2,TTCLR,TTFIN,TTRIB(25)  ,  TTSET 

CONNON  /SYSTE V  NODES,  NTI  P E,  DH  2 (  11)  ,  DH 3  (  1 1)  ,  GAIN  (  11)  ,^A0  (11), 

1  PHERR  (1 1)  ,NFD  (1 1)  ,TS  (1 1)  , CHFREQ, PHSYNP, PHSYNN , PHSY NC , N PLO HI , 

2  NNOP,TPI,N  BPO,NFSTRT, NF  STP,NFCOR,DTSAVE,DTF  ,IDECSN,NS  EF, INIT, 


*  ■-  ***-.'* -j;- 


non 


3  NNCMT,QP2,  NUMB  (1000)  ,  NSEQ  ( 1000)  ,IFILBT,  NCMT,F2,P3,  NBTJ  ( 1 000) 

C 

TP I  =  6.283185 
INCNT=0 
NCBDR= 5 
MP8NT=6 
C 

CALL  GASP 
STOP 
C 

END 

BLOCK  DATA 
BEAL*  1 6  QF2 

COMMON  /SYSTEM/  MODES,  NTYPE,  DB  2(  1 1)  ,  DB 3 (  1 1)  ,  GA IN  (1  1 )  ,TAU  (1 1 )  , 

1  PHERR  (1 1 )  ,  NPD  (1 1 )  ,TS  ( 1 1) , CHFREQ, PHSYNP, PHSYNN , PHSYNC , NPLOMI , 

2  MNOP,  TPI,  NBPO,  NPSTRT,  NFSTP, NFCUR,DTSAVB,DTP  ,  I  DECS  N,NSPF,INIT, 

3  NNCNT  ,  QP2  ,  NUNB(IOOO)  ,NSEQ(1000)  ,  N  PILBT,  MCMT,  F2,  F3,  NBTJ  { 1 000) 

DATA  NSEQ/1 000*0/, NUMB/ 1000*0/, NBTJ/ 1000*0/ 

END 

C*****SU BROOTINE  EVMTS ***********************************  *************** 

C 

C  PORPOSE 

C  PROCESS  SIMULATION  EVENTS 

C 


C 

DESCRIPTION 

OF  PARAMETERS 

C 

NSYNC 

- 

POINTER  TO  NODE  N  SYNC  PRASE  STATE. 

c 

NPHSE 

- 

POINTER  TO  NODE  N  OUTPUT  PHASE 

STATE. 

c 

NSRCE 

- 

POINTER  TO  NODE  H  SOURCE  PHASE 

STATE. 

c 

nint 

- 

POINTER  TO  NODE  N  INTERMEDIATE 

STATE. 

c 

PHSRCE 

- 

PHASE  OF  NODE  N  SOURCE. 

c 

PHSYNC 

- 

PHASE  OF  NODE  N  SYNC. 

c 

PHOUT 

- 

PHASE  OF  NODE  N  OUTPUT. 

c 

PHERR  N 

- 

PHASE  OF  NODE  N  ERROR. 

c 

PHTRK 

- 

PHASE  OF  TRACKING  LOOP  ERROR. 

c 

F  NODE 

- 

FREQUENCY  OF  NODE  N  VCO. 

c 

FSLEM 

- 

RATE  OF  CHANGE  OF  FREQUENCY  OF 

NODE  N  VCO. 

c 

FTRK 

- 

FREQUENCY  OF  TRACKING  NODE  VCO. 

c 

c 

FSLEHT 

— 

RATE  OF  CHANGE  OF  FREQUENCY  OF 

TRACKING  NODE 

V. 

c 

SUBPROGRAMS 

REQUIRED 

C  STATE, PHASE  ♦  GASP  S  OBPROGR AMS 

C 

c *************** ******************************************** ************ 

C 

SDBRODTI  NE  EVNTS(IX) 

SUBPROGRAM  DATE  780417 

DIMENSION  XX(10) 

RE  A  L*  1 6  QHFREQ,  AA1  ,  AA2,  QF2 

COMMON  /GCOH  1/  ATRIB(25)  ,JEVNT,HFA  ,HFE  (100)  ,HLE(100)  ,MSTOP,  NCRDR,  N 
1  NAPO,NNAPT,  NNATB,  NNFIL,  NNQ  (100  )  ,  NNTBY  ,  N  PRNT,  PP  ARM  (50,4)  ,TNOH,TTBES 
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2,TTCLB, TTFII, TTBIB  (25)  ,  TTSET 

COHHOH  /GCOH2/  DD  (100)  ,  DDL  (  100)  ,  DTFOL,  DTNOH ,  ISEES, L FLAG  (50)  , NFLAG , 
1  HHEQD, RMEQS , HHEQT,SS  (100)  ,  SSL ( 100) , TTHEX 
COHHOH  /GCOH 3/  AAEBB ,DTHA X, DTHIH ,DTSAV ,IITES , LLEBB, LLS AV ,LLS EV , HBE 
IBB  TTLAS  TTS AT 

COHHOH  /GCOH  4/  DTPLT  (1 0)  ,HHLOI  (25)  ,HH§ID  (25  )  ,IICBD,  IIT  AP(  10 )  ,  JJCEL 
1  (500)  ,LL  ABC  (25,2)  ,  LLABH  (25,  2) ,  LLABP  (11,2),LLABT  (25,2),LLPHI  (10)  ,LL 
2PL0  (1 0 )  ,  LLPLT,  LLSOP  (15)  ,  LLS  Y  H  (  10)  ,HHPTS,  HHCEL(25)  ,NNCLT,NNHIS,NNPL 
3T,  HNPTS(IO)  , NHSTA,NHVAB  (10)  ,PPHI  (10)  ,PPLO(10) 

COHHOH  /GCOH6/  EEMQ(IOO)  ,  IIHN(  100)  ,KKRNK(100)  , HHAZQ  (100)  ,QQTIH(100 
1)  ,SSOBV(25,5)  ,  SSTPV(25,6)  ,  VVNQ  (100) 

COHHOH  /SYSTEV  NODES,  NTY PE,  DS  2  ( 1 1)  ,DH3(11)  ,GAIH(11)  ,TAU{11)  , 

1  PHEBR  (1 1)  ,HFD  (11)  ,TS  (11)  ,CHFREQ,  pHSYNP, PHSYNN , PHSY NC,  NPLOHI, 

2  HNOP,  TP  I,  H  BPO  ,  NFSTRT,  NFSTP,NFC(JH,DTSAVE,DTF  ,1  DECS  N,RSPF,INIT, 

3  NHCNT,QF2,  HOHB(IOOO)  ,  MSEQ(IOOO)  ,  NFILBT,  NCNT,  F2,  F3,  NNTJ  ( 1 000) 

C 

C  ***** 

C  * 

C******  DETEBHIHE  HODE  FOB  EVENT  ACTIVITY  AND  ESTABLISH 

C  *  COBBECT  POINTEBS  FOB  PHASE 

C  * 

C  ***** 

C 

DATA  HHCHT/0/, NNFLAG/0/ 

N  =  ATRIB(  3) 

IF  (IX  .EQ.  4)  N  =  NODES 
CALL  STATE 

CALL  PHASE  (NT YPE,  N,  NSY  NC,  NPHSE,  NINT,  PHSYNC,  P HOOT, P  HINT, PH SRCE) 

PHSYNP=PHSYNN 

PH  SYMN  =P  HSYNC 

IF (NCNT. EQ- 1)  NPS  EQ=0 

IF(NCNT.GT.  1)  NCHTH  1=NCHT- 1 

IF  (NCNT.GT.  1)  NPSEQ=  NSEQ(NCNTH  1) 

IF ( NNFLAG. EQ . 1 -AND -I X.NE .2)  GO  TO  4012 
IF  (NMFLAG.  EQ.  1)  GO  TO  3001 

IF(NFCOH -LT.  3-OB  .  AB  S  (PHSYN  N)  -GE.ABS  (PHSYNP)  )  GO  TO  40 \2 
3001  NNF  LAG=0 

IF  (NSEQ  (NCNT)  -EQ.O)  NNFLAG=1 

IF  (HSEQ(NCNT).  EQ.O)  GO  TO  5568 

AA1  =  (QF2*1.0 Q3)  /QFLO AT ( NBPO* NSEQ  (NCNT)) 

IF  (NHTJ  (NCNT) .  EQ.O)  GOTO  310 
AA2=  1.  0562  5Q0/QFLOAT  (N  NTJ (NCNT)  ) 

GO  TO  320 
310  AA2-0.  QO 

320  QHFBEQ=  (QF  2*1 . 0Q3*QFL0 AT (NPLOHI)  *AA1*  (1 . 0Q0  ♦  AA2  )  )/1 . 0  Q3 
CHFBEQ-SNGLQ  (QHFBEQ) 

F3=CHFREQ 
GO  TO  5569 

5568  CHFHEQ  =  F2 

5569  F3*CHFBEQ 
NNODES=NODES  *  1 

C  WHITE  (NPBNT, 5045)  TNOH , CHFREQ,  F2 
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5045  FOBS AT ( •  THOM** ,E1 4.7, 51, 'CHPREQa' ,  1PB1 4.7 ,5X,' F2=* , El 4. 7) 

C  IF  (THOM.  LB. 5.  B-4)  MBIT  B(NPBHT,  4054)  NCHT,CHFREQ,MSEQ  (NCMT)  , 

C  IMHO DBS 

4054  FOB  BAT  ( •  IM  EVHTS.  ..  .  IX=4.  .  .  .HCBT=' ,  13,  21,  •  CHFREQa*  .  220.7,  21, 
1'MSBQ  ()=,,I5,/2X#»  MIODES=' ,13) 

00  14110  IJKI0B=1, HMODBS 

IF  (IJK 10 B.  LB.  MODES)  DM2  (UKMOH)  *  1.  E6*  (F2-CHFBEQ)  •Ir-I 
14110  DM3  (IJKM0B)=1.  B6*  (F3-CHIBEQ)  *TPI 
C 

C  ***** 

C  * 

c***«*«  PBOCESS  EVENT  CODE  IX 

C  ♦ 

c  ***** 

c 

C4012  IF  (TM0H.LE.5.  B-4)  MB  IT  E  (NPBNT,  4050)  CBFBEQ,  THOM 

4050  FOBS  AT  (/•  IM  ETMTS . CRFREQ  =  *,220.6,21, •  THOM  =  *,£16. 5) 

4012  GO  TO  (1000,2000,3000,4000,5000),  IX 
C 

c  ***** 

c  * 

C******  SAFE 

C  * 

C  ***** 

C 

1000  COMTIHOE 

C  IF  (TNOi.  LE.5.  B-4)  H BIT E (NPBNT, 4055)  IDECSN, P HSYMC,P HOUT,PHSRCE , 

C  1 PHERRM , P  HTBK 

4055  FORH AT ( •  IM  EVNTS. .1 X=1 * ,2 X,* IDBCSN=* ,12 ,2 X ,• PHSYNC=* , El 3. 5, 
12X,*  PHOUT=*  ,E13.  5,/,  •  PRSRCEa*  ,  E13.5,  2X,*PHEHRN  *,E13.5,2X, 

2*  PHTBK=* ,E13.5) 

C 

C******  SAVE  DATA  FOB  PLOTTING  AND  STATISTICS  COLLECTION 

C  *  AT  REGULAR  INTERVALS 

C 

TSAVE  =  TNOM 

PHSBCE  =  PHSRCE  *  360.0  /  TPI 

PHSYMC  *  PHSYNC  *  360.0  /  TPI 

PHOUT  *  PHOUT  *  360.0  /  TPI 

PHERRM  =  PHSYNC  -  PHOUT 

PHTRK  *  PHOUT  -  SS  (NNEQD-1)  *  360.  0  /  TPI 
FHODE  *  DD(NPHSE)  /  TPI 
FTRK  *  DD(MNEQD-I)  /  TPI 
C 

IF  ( IM  IT  .EQ.  0)  GO  TO  1010 


C 

C******  FIB  ST  SAVE  TIHE 

C  *  CLEAB  STATISTICAL  ABBAYS  AMD  INITIALIZE 

C 

CALL  CLEAR 

SSTPV (  1,  6)  =  PHOUT 

SSTPV  (4  ,6 )  *  PHEBBH 


SSTPV  (5 ,6)  =  PHTRK 
SSTPT  (6,6)  =  PHS8CE 
C 

1010  CONTINUE 

c 

IF  (NM  PIT  .EQ.  0}  60  TO  1100 
C 

c******  STOHE  DATA  POR  PLOTTI EG 

C 

XX  (  1)  =  1DECSN 
XX  (2)  =  PHSYNC 
XX  (3)  =  PHOUT 
XX  (4)  »  PHSBCE 
XX  (5)  =  PHERRN 
XX  (6)  =  PHTRK 
IDECSN  =  0 
NNC  NT=  NNCNT* 1 

IF(HN3P.EQ.  1)  CALL  GPLOT  (XX,  TNOI  ,1 ) 

IF  (BOD  (NNCNT, NNOP).  EQ.  1)  CALL  GPLOT (XX, TNON ,  1) 

C 

1100  CONTINUE 

IF  (NNSTA  .  GT.  0)  CALL  T  INS  A  ( PRO  OT  ,  TNOW  ,  1) 

IF  (NNSTA  -GT.  3)  CALL  TIHSA  (PHERRN,TNOi,  4) 

IF  (NNSTA  .GT.  4)  CALL  TIHSA  (PHTRK, TNOW  ,5) 

IF  (NNSTA  .GT.  5)  CALL  TIHS  A  (PHS  RCE,  TNOW,  6) 

IF  (NN  HIS  .GT.  2)  CALL  HISTO  (PHOOT, 3) 

IF  (NNHIS  .GT.  3)  CALL  HISTO  (PHERRN  ,4) 

IF  (NNHIS  .GT.  4)  CALL  HISTO  (PHTRK ,  5) 

IF  (NNHIS  .GT.  5)  CALL  HISTO  (PH SRCE ,6 ) 

C 

C******  SCHEDULE  NEXT  FREQUENCE  SAVE  TIHE  TO  GET  FS LEU 

C 

ATRIB(I)  =  TNOf  ♦  DTF 
AT  RIB  (2 )  =  5.0 
CALL  FILER  { 1) 

RETURN 

C 

C  ***** 

C  * 

c******  DCRCT 

C  * 

C  ***** 

C 

2000  CONTINUE 

IF  (N  .EQ.  NODES)  NFCUR  *  NFCUR  ♦  1 

C  IF (TNON. LE.5.E-4)  H RlT E ( NPRNT,  4056)  N , NO DES, N FCUR 

4056  FORH  AT  (•  IN  EVNTS.  .  .  11  =  2.  .  .  N=  •  ,  I  3,  2X ,  *NODES  = \  1 3,  2X  ,  'NFCUR  =  •  , 
113) 

IF  (NFCUR  .GT.  HFSTP)  HSTOP  =  -1 
C 

C******  SCHEDULE  NEXT  DBCISION  FOR  CORRECTION 

C 
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ATBIB(I)  =  TlOH  ♦  IBPO  *  TS  (N) 

CALL  FILER (1) 

C 

IDECSN  =  1 

IF  ( (R  FCUR  .LT.  HPSTHT)  .OB.  (H  .  NE .  BODES))  SO  TO  2500 
C 

c******  CBEATE  BSPF  SAFE  TIRES  FOB  THIS  PBARE 

C 

ATBIB(  2)  -  1.0 
ATBI B  (3 )  =  H 
ATB IB (  4)  =  0.0 
DO  2100  I  =  1,NSPF 

ATB  IB  (  1)  =  TBOi  ♦  (I  -  1)  *  DTSAFE  ♦  NFD  (K)  *  TS(N) 

CALL  FILEM(I) 

2100  CONTINUE 
C 

C******  TEST  FOB  1  BIT  ERROB  BETSEEN  SOUBCE  CHANNEL  ABD 

C  *  FREQUENCY  SOURCE  P  BAS  E 

C 

2500  CONTINUE 

EBB  =  PH  SYNC  -  PHSBCE 
IF(ABS  (EHR)  .LT.TPI)  RETURN 
C 

C******  SCHEDULE  PHASE  COBBECTIOB 

C  *  ATTBIBUTE  4  IS  THE  CORRECTION  TO  BE  BADE 

C  *  CORRECT  PHASE  IMMEDIATELY  IF  NFD(N)  =  0 

C 

ATRI B  (4  )  =  SIGN  (TPI,  ERR) 

IF  (NFD(N)  .EQ.  0)  GO  TO  3000 
ATRI  B  (1  )  =  TNOW  ♦  NFD(N)  *  TS(N) 

ATRIBf  2)  =  3 
CALL  FILER  (1 ) 

RETURN 

C 

c  ***** 

C  * 

c******  CRCTN 

C  * 

C  ***** 

C 

3000  CONTINUE 

C  IF  (TNOB.  LE.5.  E-4)  HRH  E  (NPRNT,  4057) 

4057  FORM  AT ('  IN  EFN TS. . . IX=3 » ) 

C 

C******  CORRECT  PHASE  OF  SYNC  FREQUENCY  FOR  NODE  N 

C 

C  IF(TNON.LE.5.B-4)  NBITE  (NPRNT,4999)  NCNT,  NUMB (  NCNT) 

4999  FORMAT ( '  NUHB(,#I3,')=,,I5) 

ANU  NBB=  NUMB  (NCNT) 

SS(NSYNC)  =  SS(NSYNC)  -  ANUHBH  *  ATB  IB  ( 4) 

NCNT=NCNT*1 

IF(NCNT.GT.NFILBT)  NCNT=1 
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BETOBN 


C 

C  ***** 

C  * 

C******  STATE  B VENT 

C  * 

C  ***** 

C 

4000  CONTINUE 

C  IP  (TN09.  LE.  5.  E-4)  8  BIT  E  (NPBNT,  405 1) 

4051  POBH  AT  {•  IN  EVNTS. .. .1 X=4  • ) 

PHOOT  =  PHODT  *  360.0  /  TPI 

IP  (LFLAG(I)  .  EQ.  0)  GO  TO  4010 
C  PHASE  HAXIHOH 

IP  (MNHIS  .GT.  0)  CALL  HISTO(PHOOT,  1) 

IP  (NNCLT  .GT.  0)  CALL  COLCT  (PHOOT,  1) 

C  IP  (TNOB.  LE.5.  E-4)  8  BITE  (HPBNT,  4052) 

4052  POBH  AT  {  *  IN  EVNTS.  .  .  .1  X=4  .  .  .  PH  A  SE  HAX  ') 

GO  TO  4020 

C 

4  010  CONTINUE 

IF  (LF LAG  (2 )  .  EQ.  0)  GO  TO  4020 
C  PHASE  H  IN  III  OH 

IF  (NNHIS  .GT.  1)  CALL  HISTO(PHOOT,  2) 

IP  (NNCLT  .GT.  1)  CALL  COLCT  ( P HO OT ,  2) 

C  IP  (TNO  8.  LE.5.  E-4)  8  BITE  ( NPBNT,  4053) 

4053  POBH  AT  (  '  IN  EVNTS. .. .IX=4. . .PHASE  BIN*) 

C 

4020  CONTINUE 
BETUBN 

C 

C  ***** 

c  * 

C******  PBEQUENCT  SAVE  EVENT 

C  * 

c  •••*• 

C 

5000  CONTINUE 

C  IP  (TNOH.  LE.5.  E-4)  8R  IT  E  (NPBNT,  4058) 

4058  POBH  AT ( '  IN  EVNTS. ,.IX=5’) 

PSLE8  =  (DD(NPHSE)  /  TPI  -  FNODE)  /  (TNOH  -  TSAVE) 
PSLE8T  =  (DD(NNEQD-I)  /  TPI  -  PTRK)  /  (TNOH  -  TSAVE) 
IF  (INIT  .EQ.  0)  GO  TO  5010 
INIT  =  0 

SSTP V (  2,  6)  =  PNODE 

SSTPV  (3  ,6 )  =  PSLE8 
SSTPV  (  7,  6)  =  PTBK 

SSTPV  (8  ,6)  =  PS  LENT 
C 

501  0  CONTINUE 

IF  (NNPLT  .  EQ.  1)  GO  TO  5100 
C 
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c******  STORK  DATA  FOR  PLOTTI MG 

C 

XX  l  1)  =  FNODE 
XX  (2)  =  PSLEB 
XX  {  3)  =  PTH  K 

XX  (4  )  =  PSLEiT 
BBCHT^HHCHT* 1 

IP  (NMOP.  EQ.  1)  CALL  GPLOT  (XX  ,  TS  A?  E,  2) 

IF (NOD (SBCBT, BNOP)  .EQ.1)  CALL  GPLOT  (XX, TSATE,2) 

C 

5100  CONTINUE 

IP  (NNSTA  .GT.  1)  CALL  T  INS  A  (PNODE,  TS  AT  E,  2) 

IP  (NNSTA  .GT.  2)  CALL  TIHSA  (P  SLER ,TSAVE,3) 

IP  (NNSTA  .GT.  6)  CALL  T  INS  A  (FTRK,TSATE,  7) 

IP  (NNSTA  .GT.  7)  CALL  TIflSA  (PSLE«T,TSATE,  8) 

IF  (NNHIS  .GT.  6)  CALL  HI STO  (P  NODE  ,7 ) 

IP  (NNHIS  .GT.  7)  CALL  HISTO  ( PSL  ES  ,  8) 

RETURN 

C 

END 

C*****SUBROUTI NE  I  NT LC ****************** *********  *********************** 

C 

C  PURPOSE 

C  INITIALIZE  SYSTEM 

C 

C  REHARKS 

C  CODE  IS  INCLUDED  TO  INITIALIZE  SYSTEM  IN  STEADY  STATE. 

C 

C  SUBPROGRAMS  REQUIRED 

C  PHASE  ♦  GASP  SUBPROGRAM 

C 

C***** ********************************************  ********************** 

C 

SUBROUTINE  IN  TLC 
C 

C  SUBROUTINE  DATE  770221 

DIMENSION  PSINC(11),POUT(11) 

REAL*1 6  QF2 ,AA1,AA2,QHFBEQ 
C 

COMMON  /GCOM1/  ATRIB(25)  ,  J  ET  NT,  MPA,  H  FE(  100) , NLE  ( 100 ) , MSTOP , NCRDR ,N 
1NAPO,NNAPT,NNATR,NNPIL,NNQ (1 00) ,  NNTR  Y  ,  N  ERNT  ,  PP  ARM  (  50,  U)  , TNOH, TTBEG 
2,TTC LR,TTFI N, TTRI B ( 2 5)  , TTSET 

COMMON  /GCOM2/  DD  (100)  ,  DDL  ( 100)  ,  DTFUL,  DTNON,  ISEES,LFLAG  (  50)  ,  NF  LAG  , 
1NNF.QD,NNEQS,NNEQT,SS  (100)  ,  SSL  ( 100)  , TTNEX 
COMMON  /GCOH  3/  AAERR ,DTHAX,DTMIN,DTSAT,IITES ,LLERR, ILS AT, LLS EV, RRE 
1 RR ,TTL AS ,TTS  AT 

COMMON  /GCOM  4/  DTPLT(IO)  ,HHLOH  (25)  ,  H  H  HI  D  (25)  ,IICRD,  IIT  AP  ( 10  )  ,  JJCEL 
1  (500)  ,  LLABC(25,2)  ,  LL  ABH  (25,  2)  ,  LL  ABP  ( 1 1,  2) ,  LL  ABT  (25,  2 )  ,  LLPHI  (10)  ,LL 
2PLO  (10)  ,LLPLT,LLSUP  (15)  ,LLSYM(  10)  ,  MMPTS,  NNCEL  (  25)  ,NScLT, NNHIS, NNPL 
3T,  NNPTS  (  10)  ,  NNSTA,  NNTAR  (10)  ,PPHI  (10)  ,PPLC  (10) 

COMMON  /SYSTEM/  NODES,  NTYPE,  DH  2 ( 1 1)  ,  DH  3  (  1 1)  ,  GAIN  ( 1 1 )  ,TAU  (1 1)  , 

1  P  HERB  (1  1)  ,NPD  (11)  ,TS  (1  1)  ,  CH  FR  EQ,  PHSY  NP  ,  PHS  YNN  ,  PHSY  NC  ,  NPLO  HI , 
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2  MNOP,TPI,NBPO,NFSTRT,  NFSTP,  NFCDR,  OTS  AT  E,  DTF,IDECSN ,NSPF,INIT, 

3  NNCNT  ,QF2  ,  NUMB  (1000),  NSEQ(IOOO)  ,  NFILBT,  NCNT,F2,  F3, N BTJ ( 1 000) 

C 

BRITE (NPRNT,9700) 

9700  FORMAT  (•  1*  ) 

C 

NCNT= 1 
PHSYNN=0. 

NNEQD  =  0 
NFCUR  =  1 
HCRDR=NCRDR-1 
LCR  DR=NCRDR- 2 

READ  (NCRDR,  101)  NFSTRT,  NFSTP,  NSPF 
10  1  FORM  AT  (  3110) 

READ  (LCBDR,*)  MNOP,  NPLOMI 
HRITE(NPRNT,9 100)  NFSTR T, NFSTP, NSPF 
9100  FORMAT  (10X,  '  DATA  COLLECTION  ST  ARTS  AT  FRAME*, 14, 

1  IX, 'AND  STOPS  AT  FRAME  ',15, IX, 

2  'KITH  ',14,'  SAMPLES  PER  FRAME  ') 

C 

C  ***** 

c  * 

c******  READ  INPUT  PARAMETERS  AND  CONVERT  TO  SIMULATION 

C  *  INPUT  FORMS 

C  * 

C  ***** 

c 

C  ******  ENTER  CHANNEL  FREQUENCY  IN  MHZ  AND  NUN  BEP 

C  *  OF  BITS  per  OPPORTUNITY 

C 

READ  (NCRDR, 102)  CHFREQ, NBPO 

102  FORMAT  (F  10. 0,110) 

C 

C* ** ***  ENTER  NUMBER  OF  NODES  AND  TYPE  OF  SYSTEM 

C 

REA  D  (NCR  DR,  103)  NODES, NTYPE 

103  FORMAT  (2110) 

WRITE(NPRNT, 9120)  NODES, NTYPE 

9120  FORMAT  (///,1  OX, 'SYSTEM  PARAMETERS  FOR  •  ,  1 3,  IX  ,  •  NODES'  / 

1  10X , '  ALL  NODES  ARE  TYPE', 12,1  X, 'SYSTEMS') 

RE  AD (MCR DR, 110)  NFILBT 
110  FORMAT  (110) 

READ  (MCRDR,  120)  (NSEQ(IAD)  ,IAD=1  ,NFIIBT) 

READ  (MCRDR,  120)  (NWTJ(IAD)  ,  I  AD=  1,  N  FIL BT) 

REA  D (M  CR  DR,  120)  (NUMB (IAD)  ,IAD  =  1 , NFILBT) 

120  FORMAT  (1  415  ) 

BR  IT  E  (  6,  140) 

140  FORMAT  (/10X,' NS  EQ  ARRAY  LISTED') 

BRITE(NPRNT,  125)  (NSEQ  (JAL)  ,  JA  L=1  ,  NFI  LBT) 

125  FORMAT  (  10X,  1415) 

BR  ITE  (6 , 14  5  ) 

145  FORMAT  {/10X, 'NBTJ  ARRAY  LISTED') 
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WRITE(NPRNT,125)  (NWTJ  (JAL)  ,  JAL=  1,  NFILBT) 

WRITE  (6,  130) 

130  FORHAT  (/10X,*  NUHE  ARRAY  LISTED*) 

WRITE  (NPRNT, 125)  (NUHB  (JAL)  #  JAL=  1,  NFILBT) 

C 

c ******  ENTER  PARABETER  SETS  FOR  EACH  NODE 

C 

DO  1000  N  =  1  ,  NODES 
WRITE(NPBNT,9130)  N 
9130  FORHAT  (//,10X, ‘NODE*  ,  13) 

C 

c******  ENTER  DATA  TRANSFER  FR  EQDENCY  IN  HHZ,  FIXED  DELAY 

C  *  IN  BITS  AND  VCO  NOBIN  AL  FREQUENCY  IN  H  HZ 

C 

REA.,  (NCRDR,104)  QF2, NFD (N) , F3 

104  FORH  AT  (Q  10- 3, 110, FI  0-0) 

F2  =  SNGLQ  (QF2) 

IF  ( N  .N  E.  1)  GO  TO  112 

AA  1  =  (QF2*1. 0Q3)  /QFLO  AT  (NBPO*NS  EQ(1)  ) 

IF(NWTJ(  1)  .EQ.O)  GO  TO  310 
AA2=1-05625QO/QFLOAT  (NWTJ  (1)  ) 

GO  TO  320 
310  AA  2=0. 0Q0 

32  0  QH  F  BE  Q=  (QF2*  1 . 0  Q3*  Q  FLO  AT  (N  PL  OR  I)  *AA1*  (  1.  OQO+AA  2)  )  /1.0Q3 
CH  FR  EQ  =SN  GLQ (QHPREQ) 

112  F3  =CHF  REQ 

IF (N.N  E.  1)  GO  TO  111 
WRITE  (NPRNT  ,9  110)  CHFREQ,  NBPO,  HNOP 
9110  FO  RH  AT ( ///,  10  X  ,  '  T 1  FREQUENCY  SOURCE  = • , 1  PEI  5 . 7 , 1 X , • MHZ  •/ 

1  1 0  X , ' NU  HBER  OF  BITS  PER  OPPORTUNITY  =',I4/10X, 

2  •  EVERY*  ,16,  'TH  DATA  POINT  IS  PLOTTED') 

111  WRITE(N?RNT,9140)  F2,NFD(N),F3 

9140  FORHAT  (/,10X,  'T2  FREQUENCY  SOURCE  =  •  ,  IP  E  1 5.  1,  IX  ,  •  B  HZ  '  / 

1  10X,' CORRECTION  APPLIED  IN  BIT 13, IX, 'OF  FRARE'/ 

1  10X,'VCO  NOMINAL  FREQUENCY  =  '  , 1  PE  1 5 . 7 , 1 X , ' HHZ • ) 

TS  (N)  =  1.0  /  (F2  *  1.0E6) 

DW  2  (N )  =  1.0E6  *  (F2  -  CHFREQ)  *  TPI 
DW  3  ( N)  =  1.0E6  *  ( F3  -  CHFREQ)  *  TPI 
C 

C******  ENTER  SYNC  PHASE  AND  OUTPUT  PHASE  IN  DEG 

C 

READ(N  CRDF,  105)  PSYNC  (N)  ,POUT  (N) 

105  FORBAT  (2F10.0) 

W  R  IT  E  (  N  P  RN  T  ,  9  1  4  5)  PSYNC  (N)  ,POUT(N) 

9145  FORHAT  (/,10X,'T2  INITIAL  PHASE  =  '  ,  F5. 0,  IX  ,  •  DEG'  / 

1  10X,'NODE  OUTPUT  PHASE  =  '  ,F5 . 0 , 1  X  ,  'DEG  •) 


PSYNC.N)  =  PSYNC(N)  *  TPI  /  360.0 
POUT(N)  =  POUT(N)  *  TPI  /  360.0 

C 

c******  LOOP  BANDWIDTH  IN  HZ  AND  DAMPING  (SECOND  ORDER  ONLY) 

C 


IF  (NTTPE  .EQ.  2)  GO  TO  100 


E-ll 


o  n  o  on^nnnnnnnnn  n  non  n  -»  o  o  n 


FIRST  ORDER  LOOP 

READ  (NCRDR, 106)  FN 
06  FORM  AT  (F  10.0) 

WRITE  (NPRHT,9150)  FN 

9150  FORMAT (/, 10X,' NATURAL  FREQUENCY  (FN)  =  •  ,1  PE15.  7 , 1 X  ,  •  HZ  • ) 

WN  =  T  PI  *  FN 

CONVERT  BANDWIDTH  INPUT  TO  GAIN  FOR  SIMULATION 
GAIN  (N)  =  WN 
NNEQD  =  NNEQD  ♦  2 
GO  TO  200 

SECOND  ORDER  LOOP 

100  CONTINUE 

READ  {NCRDR, 105)  FN,Z  ET  A 
HRITE(NPRNT,9  160)  FN  ,ZETA 

9160  FORMAT  (/,1  0  X,  '  NATURAL  FREQUENCY  (FN)  =•  ,  IP  E  1 5.  7,  IX  ,  •  HZ  •  / 

1  10X, 'DAMPING  (ZETA)  =  '  ,  0  PF  j.  3) 

WN  =  TP  I  *  FN 

CONVERT  BANDWIDTH  INPUT  AND  DAMPING  TO  GAIN  AND  TAU  FOR  SIMULATION 

GAIN  (N)  =  WN  *  WN 

TAU  (N)  =  2  *  ZETA  /  WN 

NNEQD  =  NNEQD  +  3 

200  CONTINUE 

1000  CONTINUE 

***** 

* 

******  TRACKING  LOOP  PARAMETERS 

* 

***** 


******  ENTER  VCO  NOMINAL  FREQUENCY  IN  MHZ 

READ(NCRDR,  106)  F3 
F3  =CHF  HEQ 

WRITE(NPRNT,9170)  F3 

9170  FORMAT  (///,10X, 'TRACKING  LOOP  PARAMETERS'/ 

1  10X,'VCO  NOMINAL  FREQUENCY  = »  ,  1  PE  1  5 . 7 , 1  X,  '  MHZ ’) 

DW3  ( N)  =  1.0E6  *  ( F3  -  CHFREQ)  *  TPI 

******  ENTER  LOOP  BANDWIDTH  IN  HZ  AND  DAMPING 

READ  (NCRDR, 105)  FN,Z  ET  A 
WRITEJNPRNT,  91  60)  FN,ZETA 
WN  =  TPI  *  FN 

C  CONVERT  BANDWIDTH  INPUT  AND  DAMPING  TO  GAIN  AND  TAll  FOR  SIMULATION 
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nnn  n  non  nonn  nnnnono 


GAIN  {H )  =  IN  *  IN 
TAO(N)  =  2  *  ZETA  /  BN 
NNEQD  =  NNEQD  ♦  2 
NN  EQT  =  NN  EQD 

***** 

* 

******  INITIALIZE  SYSTEM  ETENTS 

* 

***** 


DO  2000  N  =  1 ,  NODE  S 

CALL  PHASE  {NTY  PE,  N,  NSY  NC,  N  PHS  E,  N  INT,  P  HSYNC,  P  BOUT,  P  HINT, PH  SRCE) 
IF  (NTIPE  .EQ.  2)  GO  To  1100 

FIEST  ORDER  SYSTEM 


INITIALIZE  PHASES 
SS(NSYNC)  =  PHSRCE  *  PSYNC(N) 
SS(NPHSE)  =  SS(NSYNC)  ♦  POOT(N) 
GO  TO  1200 


SECOND  ORDER  SYSTEM 


1100  CONTINUE 

INITIALIZE  PHASES 

SS  { N ST N C)  =  PHSBCE  ♦  PSYNC(N) 

SS(NPHSE)  =  SS(NSYNC)  ♦  POUT  (N) 

SS  (N INT )  =  -DB3  (N)  /  GAIN(N) 

1200  CONTINUE 

PHERR(N)  =  0.0 

******  SCHEDULE  FIRST  DCRCTN  TIME 


AT  RIB  {1  )  =  0.0 
ATR  IB  (  2)  =2.0 
ATRIB  (3 )  =  N 
AT R I B (  4)  =  0.0 
CALL  FILEM  (1) 

2000  CONTINUE 
C 

N  =  NODES  ♦  1 
NINPUT  =  NPHSE 
NPHSE  =  NNEQD  -  1 
NINT  =  NNEQD 
SS (NPHSE)  =  SS  (NINPUT) 

SS  (N  INT)  =  -DH3  (N)  /  GAIN  (  N) 

IDECSN  =  0 
I  NIT  =  1 

DTSAVE  =  FLOAT  (NBPO)  *  TS(NODES)  /  FLOAT  (  NS  PF) 
DTPLT(I)  =  DTSAVE*FLCAT  (NNOP) 

DTPLT  (2  )  =  DTSAVE*FLOAT  (HBOP) 
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DTP  =  0.1  *  TS  (NODES) 

RETURN 

END 

C*****SUB ROUTINE  PHASE*** **************************************** ******* 

C 

C  PURPOSE 

C  DETERMINE  POINTERS  AND  CURRENT  VALUES  POR  NODE  STNC  AND 

C  OUTPUT  PHASE  BASED  ON  TYPE  OF  NODE  AND  NODE  NUMBER 

C 

C  SUBPROGRAMS  REQUIRED 

C  NONE 

C 

c *****************  ******************  ********* ******************** ******* 

c 

SUBROUTINE  PH  AS  E  (NTY  PE,  NODE,  NS  Y  NC,  NPESE, NINT,P HSYNC , PHOUT, PHINT, 

1  PHSRCE) 

C 

C  SUBROUTINE  DATE  761214 

C 

COMMON  /GCOM  2/  DD(100)  ,DDL(1QO)  ,DTF UL ,DTNOi ,IS EES , LFLAG (50 )  ,  NFL  AG, 
1 NNE QD, NNEQS , NNEQT, SS (100) , SSL { 100) , TTN EX 
C 

PHSRCE  =  0.0 
C 

GO  TO  (100,200)  ,  NT  Y  PE 
C 

C  FIRST  ORDER  SYSTEM 

C 

100  CONTINUE 

NSYNC  =  2  *  NODE  -  1 
NPHSE  =  NSYNC  ♦  1 

NS  BCE  =  NSYNC  -  1 

P HSYNC  =  SS (N  SYNC) 

PHOUT  =  SS  (NPHSE) 

IF  (NSRCE  .GT.  0)  PHSRCE  =  SS(NSRCE) 

RETURN 

C 

C  SECOND  ORDER  SYSTEM 

C 

200  CONTINUE 

NSYNC  =  3  *  NODE  -  2 
NPHSE  =  NSYNC  ♦  1 

NSRCE  =  NSYNC  -  1 

NINT  =  NPHSE  ♦  1 
PH  SYNC  =  SS  (NSYNC) 

PHOUT  =  SS(NPHSE) 

PHI  NT  =  SS(NINT) 

IF  (NSRCE  .GT.  0)  PHSRCE  =  SS(NSRCE) 

RETURN 

C 

END 

C  *****  SUB ROD TI NE  SCOND**********  **************************************** 
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C  PURPOSE 

C  DETECT  MINIMUM  ADD  HAIIMUH  PHASE  PEAKS  AT  LAST  NODE  IN  CASCADE 

C 

C  BE  HARKS 

C  PHASE  MIN  AND  HA X  DETECTED  SITHIN  TNG  DEGBEES. 

C 

C  SDBPROGBAHS  REQUIRED 

C  GASP  SUBPROGRAMS 

C*************************************************  ********************** 

C 

SUBROUTINE  SCOND 
SUBROUTINE  DATE  770110 
REAL  *  1 6  QF2 

COHHON  /GCOH1/  ATRIB(25)  ,  J  EV  NT,  HFA,  MFE(  100)  ,  ML  E  ( 1 00)  ,  MSTOP  ,  NC  RDR  ,  N 
1NAPO,NNAPT,NNATR,NNFIL,NNQ  (100)  , NNTB Y , N PENT , PP ARM ( 5 0, 4)  f T  NOB ,  TTB  EG 
2,TTCLB,TTFI  N,TTRIB(2  5)  ,  TTSET 

COHHON  /GCOH2/  DD(100)  ,DDL(100)  ,  DTFUL,  DTNOW ,  IS  EES,  L  FLA  G  ( 50)  ,NFLAG, 
1NNEQD,NNEQS,NNEQT,SS  (100),  SSL  (100)  , TTNEX 
COHHON  /SYSTEM/  NODE  S,NTYPE,DW2  (1 1  )  ,DH3  (1 1 )  ,GAI  N  (1 1  )  ,T  AU  (1  1)  , 

1  PHERR  (11)  ,NFD(11)  ,TS(11)  ,  CHFR EQ,  P HSYNP, P HS INN , PHSYNC, NPLOHI , 

2  HNOP,  TP  I,  NBPO,NFSTRT,NFSTP,  NFCUR  ,DTSA  VE,DTF  ,1  D  ECS  N,  NS  PF,  IN  IT, 

3  NNCNT,QF2 , NUHB(IOOO) , NS EQ (1 000 ) , N FILBT, NCNT, F2, P3, NWTJ (1000) 

GO  TO  (100  ,200)  ,  NT Y  PE 

100  CONTINUE 

FIRST  ORDER  SYSTEM 
NPHSE  -  l  *  NODES 
GO  TO  300 
200  CONTINUE 

SECOND  ORDER  SYSTEH 
NPHSE  =  3  *  NODES  -  1 
300  CONTINUE 

LF  LAG  (1  )  =  KROSS  (-NPHSE, 0,0, 0,-1,  100.0) 

LFL AG(  2)  =  KROSS(-NPHSE,0,0,0,1 ,100.0) 

RETURN 
END 

♦  ♦♦♦♦SUBROUTINE  STATE** *************************************  *********** 

PURPOSE 

CALCULATE  SYSTEH  STATE  DERIVATIVES 
BE  HARKS 

THE  SYSTEH  IS  A  CASCADE  OF  NODES  FOLLOWED  BY  A  TRACKING  LOOP. 
ALL  NODES  IN  CASCADE  ARE  THE  SAME  TYPE  (FIRST  OP  SECOND  ORDER). 
THE  TRACKING  LOOP  IS  A  SECOND  ORDEB  LCOP. 


SUBPROGRAMS  REQUIRED 


oono  nnonnnnnnoo  n  r>  n  n  o  r>  n 


SOME 

*********************************************************************** 


SUBROUTINE  STATE 
SOBBOOTISE  DATE  770221 


BEAL* 1 6  QF2 

COHHOW  /GCOM1/  ATRIB  (25)  ,  J  EV  NT  ,  HFA,  M  FE(  100)  , ML E  ( 100 ) , MSTOP , NCRDR ,N 
1HAPO, NN APT, NNATR,NNFIL,NNQ(100)  , NHTBY, NPBNT , PPARM ( 50, 4)  ,T  MOW , TTBEG 
2,TTC  LR,TTFI N,TT  BI B (2  5)  , TTSET 

COMMON  /GCOH2/  DD(100)  ,DDL(100)  ,  DTFUL  ,  DTNOB ,  ISEES,  LFLAG  (  50)  ,  NF  LAG, 
1HNEQD,NMEQS,NNEQT,SS  (100)  ,  SSL  (  100)  , TTNEX 
COMMON  /SYSTEM/  NODES,  NTYPE,DS2  (11)  ,DN3  (11 )  , GAIN  (11  )  ,TAU  (11)  , 

1  PHEBB  (1 1 )  ,  NFD  (11)  ,TS  (11)  , CHFH EQ, P HSY NP, P HS YNN , P HSYNC, NPLOMI , 

2  MNOP,  TPI,  NBPO,  NFSTRT,  NFSTP,NFCUR,DTSAVE,DTF  ,IDECSN,NS  PF,INIT, 

3  N NC NT  , QF2  ,  NUMB  (1000)  ,  NS  EQ  ( 1000)  ,  N  F ILBT,  NCNT,  F2,  F3,  NHTJ  (1000) 


***** 

* 

******  SOLVE  FIRST  OR  SECOND  ORDER  SYSTEM  EQUATIONS 

*  FOR  ALL  NODES 

* 

***** 

IF  (NT YPE  .  EQ.  2)  GO  TO  1000 
******  FIRST  ORDER  SYSTEM  EQUATIONS 


DO  100  H  =  1  , NODES 

NSYNC  =  2*N  -  1 

NPHSE  -  NSYNC  ♦  1 

PHERB(N)  =  SS (N  SYNC)  -  SS(NPHSE) 

DD (NSYNC)  =  DH2(N) 

DD  (NPHSE)  =  GAIN  (N)  *  PHERB(N)  ♦  DM3(N) 
100  CONTINUE 


GO  TO  2000 

******  SECOND  ORDER  SYSTEM  EQUATIONS 

1000  CONTINUE 

DO  1100  N  =  1 , NODES 
NSYNC  =  3  *  N  -  2 
NPHSE  =  NSYNC  ♦  1 
NINT  =  NPHSE  ♦  1 

PHEBR(N)  =  SS  (NSYNC)  -  SS(NPHSE) 

DD (NSYNC)  =  DW  2 (N) 

DD  (NPHSE)  =  GAIN(N)  *  (TAU(N)  *  PHERR(N)  ♦  SS(NINT))  ♦  DW3  (N) 
DD  (NINT)  =  PHERB(N) 

1100  CONTINUE 


r>oo  n  n  n  r>  c>  n  n 


V  i 


***** 

* 

****** 

* 

***** 


TRACKING  LOOP 


2000  CONTINUE 


****** 


SET  OP  POINTERS  FOR  TRACKING  LOOP 


N  *  NODES  ♦  1 
NINPUT  »  NPHSE 
NPHSE  =  NNEQD  -  1 
NINT  =  NNEQD 

PHERR(N)  =  SS  (NINPUT)  -  SS  (NPHSE) 

DD  (NPHSE)  «  G  AIN  (N)  *  (TAU(N)  *  PHERR(N) 
DD  (NINT)  =  PHEBB(N) 

RETURN 

END 


SS(NINT)  )  *■  DW3  (N) 
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APPENDIX  F 


HYLBUR  PROGRAM  SOURCE  LISTING 
A.  DATA  SET  MEMBER  • R3893. LIB (NIHGO ) •  LISTING 

1.  ;  NIHGO 

2.  SET  EXE  NOL  TER 

3.  CLR  ACT 

4.  SET  *  LAST 

5.  SET  ESC  X 

6.  J  N 1  =  0 

7.  READ  STRING  SO  PROMPT  *3  CHAR  JOB  IDENT?  ' 

8.  READ  VALUE  NO  PROMPT* RUN  TIHE(SEC)?  • 

9.  *+1  //R3893XS0  JOB  ( 1  135200G,  2  GQ5,  %N0,  8,  1 92  )  ,  *  JOHN  J  CORMACK', 

10.  **1  //  NOTIFY=R3893, MSGCLaSS=Q 

11.  READ  STRING  SI  PROMPT*  DATA  DSN  MEMBER  NAME?  • 

12.  IF(*XS1*  EQ  *•)  X  39.000 

13.  a  N1  =  N 1  ♦  1 

14.  **1  //PRINTXN1A  EXEC  PGM= UPRINT, REGION=60K 

15.  **1  //SYSPRINT  DD  SYSOUT=A 

16.  *  +  1  //SYSUT1  DD  DSN=  R389  3.  DATA  (XS  1)  ,  DISP=SHR 

17.  *+1  //PROGRUN  EXEC  PGH=GSC24,  REG ION=  192K 

18.  *+1  //STEPLIB  DD  DSN=R389  3.  LOAD,  DISP=SHR 

19.  *+1  //FT 04 PO0 1  DD  SYSOUT=A 

20.  **1  //FT05F00 1  DD  DS N=  R3893.  DATA  (%S  1)  ,  DISP=S HR  , LA BEL=  (,  ,  ,  IN) 

21.  READ  STRING  S2  PROMPT* GSC24.F8XX .DATA  NUMBER?  • 

22.  IF  (»%S2  '  NE  •’)  a  S2=LPAD(S2,  2,*  0*  ) 

23.  **1  //FT06F00 1  DD  STSOUT=A 

24.  IF  ( *  %S2  *  NE  ")  **1  //  FT  08  POO  1  DD  DSN  =  R3893.  GSC24.  F8XS2.  DATA,DISP=SH 

25.  IF  (•  XS2*  EQ  *•)  *  +  1  //FT08F001  DD  SYSOUT=A 

26.  **1  //FT09F001  DD  DS  N=6TE MP1 , DIS  P=  (N  EW,  DEL  ETE)  ,  SPAC2=  (CYL,  1)  , 

27.  *>1  //  DCB=(RECFM  =  VS)  ,UNIT=SYSDA 

28.  **1  //FT10F001  DD  DS  N=6T EMP2,  DIS P*  (N  EN ,  DELETE) ,  SPACE=  (CYL  ,  1)  , 

29.  **1  //  DCB=  (RECFM=VS)  ,  UNIT=SYSDA 

3  0.  *♦  1  //FT  1  1 FO 0  1  DD  DS  N=5TEMP3,  DISP=  (N  EN,  DELETE)  ,S  PACE=  (CYL,  1)  , 

31.  *+1  //  DCB= (RECFH=V3) , UNIT=SYS DA 

32.  **1  //FT12F001  DD  DS  N=fcT  EMP4,  DIS  P=  (N  EH  ,  DEL  ET  E)  ,  SPACE=  (CYL,  1)  , 

33.  **1  //  DCB=(RECFH  =  VS)  , ONIT=SYSDA 

34.  IF  ( 'XS2 '  EQ  ••)  X  11.000 

35.  **1  //PPINTXN1B  EXEC  PGH  =  UPRINT,REGION=60  K 

36.  **1  //SYSPRINT  DD  SYSOUT=A 

37.  **1  //SYSUT1  DD  DS N= R3893. GSC24.  F8XS 2. DATA, DISP  =  SHR 

38.  X  11.000 

39.  **1  // 

4  0.  LIST 

41.  READ  STRING'SO  PROHPT'RUN  JOB(Y  OR  N) ?  ' 

42.  IF  (*XS0*  NE  * Y » )  X  46.000 

43.  FRA 

44.  SUB  UNN 

45.  JO 

46.  SET  ESC 
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B.  DATA  SET  MEMBER  *  B3893.  LIB  (GAS PH OD)  '  LISTING 


1.  ;  GASPNOD 

2.  SET  EXE  MOL  TEE 

3.  SET  ESC  6 

4.  a  SO  -  DATE 

5.  HEAD  STRING  S3  PROMPT*  DSN?  • 

6.  IP( *6  S3*  EQ  'STOP')  X  172.000 

7.  IP  (  *SS3  *  EQ  ••)  X  9.  000 

8.  B  S6  =  S3 

9.  BEAD  STBING  S4  PHOHPT'HEMBER  NAME?  1 

10.  IF  (•  6S4  *  EQ  ••)  X  62.000 

11.  OSE  6S6  (&S  4)  CLB 

12.  3  SI  =  SOBSTB  (SO,  10, 2) 

13.  d  S2  =  SOBSTB  (SO, 13,2) 

14.  a  SS  =  SOBSTB  (SO, 16, 2) 

15.  LIST  FIRST  ONN 

16.  CH  24/25  TO  *6S1*  NOL 

17.  CH  29/30  TO  *£S2*  NOL 

18.  CH  34/35  TO  *SS5'  NOL 

19.  LIST  *  ONN 

20.  CONH 

21.  LIST  15  ONN 

22.  READ  STRING  SO  PROMPT '  LOB  EB  VALOE?  • 

23.  IF  ('SSO*  EQ  •*)  X  32.000 

24.  a  SO  =  LPAD(  SO,  10) 

25.  CH  31/40  TO  *  6S0 ’  NOL 

26.  READ  STBING  SO  PROMPT* OPPER  VALOE?  • 

27.  IP  ( '6S0 '  EQ  «•)  X  30.000 

28.  a  SO  =  LPAD (SO, 10) 

29.  CH  4  1/50  TO  »6S0»  NOL 

30.  LIST  *  ONN 

31.  COMM 

32.  LIST  16  ONN 

33.  READ  STRING  SO  PROMPT*  LONER  fALOE?  * 

34.  IF  ('  GSO*  EQ  ••)  X  43.000 

35.  a  SO  =  LPAD(S0,  10) 

36.  CH  31/40  TO  *  SSO '  NOL 

37.  BEAD  STRING  SO  PROMPT*  OPPER  VALOE?  • 

38.  IF  ( ' 6S0  *  EQ  ")  X  41.000 

39.  a  SO  =  LPAD  (SO,  10) 

40.  CH  4  1/50  TO  *6S  0 '  NOL 

41.  LIST  *  ONN 

42.  COMM 

43.  LIST  22  ONN 

44.  READ  STRING  SO  PROMPT '  ABSOLOT  E  LOCAL  TRONCATION  ERROR 

45.  IFC&SO*  EQ  ")  X  48.000 

46.  a  SO  =  LPAD(  SO,  10) 

47.  CH  11/20  TO  • 6S0 •  NOL 

48.  BEAD  STRING  SO  PROMPT* RELATI VE  ERROR?  • 

49.  IF  ('  SSO  *  EQ  «•)  X  52.000 
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50.  9  SO  =  LPAD( SO, 10) 

51.  CH  21/30  TO  *$S0*  NOL 

52.  BEAD  STB IH G  SO  PBOHPT' HIH  STEP  SIZE?  • 

53.  IF('6S0*  EQ  •')  X  56.000 

54.  3  SO  =  LPAD  (SO,  10) 

55.  CH  31/40  TO  'ftSO*  HOL 

56.  BBAD  STBING  SO  PBOHPT*  HAX  STEP  SIZE?  • 

57.  IFfSSO*  EQ  •')  X  60.000 

58.  3  SO  =  LPAD (SO , 10 ) 

59.  CH  41/50  TO  *SSO*  HOL 

60.  LIST  *  OHN 

61.  COBH 

62.  LIST  25  DNH 

63.  BEAD  STBING  SO  PBOHPT • FBAHE  STABT  ROBBER?  • 

64.  IF  (*  6  SO'  EQ  ••)  X  73.000 

65.  BEAD  STBING  SI  PBOHPT' FBA HE  STOP  NOHBEB?  • 

66.  IF  (*6S1  •  EQ  ")  X  69.000 

67.  i  SI  =  LPAD  (SI,  10) 

68.  CH  11/20  TO  'ESI'  NOL 

69.  3  SO  =  LPAD  (SO, 10) 

70.  CH  1/10  TO  *SSO*  NOL 

71.  LIST  *  ONN 

72.  COBH 

73.  LIST  26  ONN 

74.  BEAD  STBING  SO  PRCNPT'NBPO?  » 

75.  IF  ('  SSO*  EQ  ")  X  79.000 

76.  a  SO  =  LPAD(SO.IO) 

77.  CH  11/20  TO  'SSO'  ONN 

78.  COBH 

79.  LIST  27  ONN 

80.  read  string  so  osing  27  col  i/io 

81.  a  SO  =  LTRIH  (SO ) 

82.  3  NO  =  SO 

83.  READ  STRI NG  SO  PROMPT' NODE  TYPE?  • 

84.  IF  ('  SSO*  EQ  ")  X  83.000 

85.  a  S7  'O' 

86.  IF  ('6  SO*  EQ  •  2*  )  a  S7  *1* 

87.  3  SO  =  LPAD(  SO,  1 0) 

88.  CH  11/20  TO  ' 6S0 '  ONN 

89.  LIST  26  ONN 

90.  IF  (NO  EQ  1)  3  S9  *28* 

91.  IF  (NO  EQ  2)  d  S9  *28,31* 

92.  IF  (NO  EQ  3)  a  S9  '28,31,  34* 

93.  IF  (NO  EQ  4)  a  S9  *  28,31,34,37* 

94.  IF  (NO  EQ  5)  3  S9  *28,31,34,3  7,40' 

95.  IF  (NO  EQ  6)  3  S9  *28,31,34,37,40,  43* 

96.  IF  (NO  EQ  7)  3  S9  *  28,31,34,37,40,43,46' 

97.  IF  (NO  EQ  8)  3  S9  *28,31,34,37,40,4  3,  46,49* 

98.  LIST  S S9  ONN 

99.  a  HO  =  28  ♦  (NO  *  3) 

100.  LIST  tV  0  ONN 

101.  BEAD  STRING  SO  PRCHPT'TI  RATE?  • 
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10  2. 

IF(  *SS0  • 

EQ 

•  • 

) 

X  110.000 

103. 

a  so 

=  LPAD (SO 

,10) 

104. 

CH  1/10  TO  • 

6S0* 

IM  26  MOL 

10S. 

CH  21/30 

TO 

•SSO 

•  IM  SS9  MOL 

106. 

CH  1/10  TO  • 

6S0* 

IM  SMO  MOL 

107. 

LIST  ; 

26  UNN 

1  08. 

LIST 

6S9 

OHM 

109. 

LIST 

SHO 

ONN 

110. 

COHN 

111. 

IF(HO 

EQ 

D 

a 

Ml 

=  28 

112. 

IF  (MO 

EQ 

2) 

a 

HI 

=  31 

113. 

IF  (HO 

EQ 

3) 

a 

Hi 

=  34 

114. 

IF  (MO 

EQ 

4) 

a 

HI 

=  37 

115. 

IF  (MO 

EQ 

5) 

a 

HI 

=  40 

116. 

IF(HO 

EQ 

6) 

a 

HI 

=  43 

117. 

IF  (MO 

EQ 

7) 

a 

HI 

=  46 

118. 

IF  (NO 

EQ 

8) 

a 

Ml 

=  49 

119. 

a  m2 

=  28 

120. 

LIST 

6H2 

UNN 

121. 

read 

STRING 

SO 

ERCH  PT  *T  2  RATE? 

i  t 

122. 

IF(  *6S0  • 

EQ 

t  • 

) 

X  125.000 

123. 

a  so  : 

=  LPAD (SO 

,10) 

124. 

CH  1/10  TO  • 

6  SO* 

UNN  IN  SH2 

125. 

IF  (M2 

GE 

Ml) 

X 

:  1 

29.000 

126. 

a  M2  : 

=  M2  ♦ 

3. 

127. 

COHH 

128. 

X  120 

.000 

129. 

COHN 

130. 

IF  (MO 

EQ 

1) 

a 

S8 

•29* 

131. 

IF(M0 

EQ 

2) 

a 

S8 

•  29,32* 

132. 

IF  (MO 

EQ 

3) 

a 

S  8 

•29,32,35* 

133. 

IF  (MO 

EQ 

4) 

a 

S  8 

•  29,32,35, 38* 

134. 

IF(NO 

EQ 

5) 

a 

S8 

•29,32,35,38, 

41 

t 

135. 

IF  (NO 

EQ 

6) 

a 

S8 

•29,32,35,  38, 

4  1 

,  44' 

1 

136. 

IF(NO 

EQ 

7) 

a 

S8 

•  29,32,35,38, 

41 

#44  ( 

r4  7  ' 

137. 

IF  ( NO 

EQ 

8) 

a 

S  8 

•29,32,35,  38, 

41 

#44, 

,47,50* 

138. 

LIST 

SS8 

UNN 

139. 

read 

STRJ 

:ng 

SO 

PROHPT*  INITIAL 

PHASE 

OF  OUTPUT? 

140. 

IF  (’ 6 

SO* 

EQ 

« • 

)  : 

X  143.000 

141. 

a  so 

=  LPAD ( 

so 

>,H 

0) 

142. 

CH  11, 

/20 

TO 

» SSO1 

'  UMM  IN  SS8 

143. 

COHH 

144. 

IF(NO 

EQ 

D 

a 

S5 

•  30* 

145. 

IF  (MO 

EQ 

2) 

a 

S5 

•30,33* 

146. 

IF(MO 

EQ 

3) 

a 

S5 

•  30,33,36* 

147. 

IF  (HO 

EQ 

4) 

a 

S  5 

•30,33,36,39' 

148. 

IF  (MO 

EQ 

5) 

a 

S5 

•30,33,36,39, 

42 

9 

149. 

IF(NO 

EQ 

6) 

a 

S5 

•30,33,36,39, 

42 

#45' 

i 

150. 

IF  (HO 

EQ 

7) 

a 

S5 

•30,33,36,39, 

42 

#45, 

48' 

151. 

IF(NO 

EQ 

8) 

a 

S5 

•  30,33,36,39, 

42 

#45 , 

r*»  8  ,5  1  * 

152- 

LIST 

SS5 

UMN 

153. 

READ 

STRI 

:ng 

SO 

PROMPT ' FM?  • 
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154.  IFfSSO*  EQ  ••)  X  158.000 

155.  d  SO  =  LPAD  (SO, 10) 

156.  CH  1/10  TO  '6S0*  IB  £S5  HOL 

157.  LIST  6S5  COL  1/10  0MB 

158.  IF  (*  6S7  '  ME  •  1  •  )  X  163.000 

159.  BEAD  STBING  SO  PROMPT* ZETA?  * 

160.  IF  (*  6SQ  *  BQ  ••)  X  163.000 

161.  »  SO  =  LPAD  (SO,  10) 

162.  CH  11/20  TO  •  BS  0 '  IN  6S5  NOL 

163.  LIST  6S5  OBN 

164.  SH  NAME 

165.  COMM  HONE  =  HO  SAVE,  STOP  =  PROGRAM  END 

166.  READ  STRING  SO  PROMPT  *  MEMBER  NAME?  * 

167.  I P ( *  6  SO  *  EQ  'NONE*)  X  4.000 

168.  IF  (  *6S0  *  EQ  ’STOP*)  X  172.000 

169.  IF  (*  6S0*  EQ  *•)  S  *  CAR  REP 

170.  IFfSSO*  NE  ••)  S  *(&S0)  CAR  REP 

171.  X  4.000 

172.  SET  ESC 

C.  DATA  SET  MEMBER  » R3893 .LIB ( GASPGO) *  LISTING 

1.  ;  GASPGO 

2.  SET  EXE  NOL  TER 

3.  SET  ESC  X 

4.  READ  STRING  S2  PROHPT*3  CHAR  IDENTIFIER?  • 

5.  I F  ( *XS2  *  EQ  ••)  X  71.000 

6.  CLR  ACT 

7.  SET  *  LAST 

8.  a  HO  =  LAST  ♦  1 

9.  READ  VALUE  N1  PROMPT*  BON  TIME  (SEC)?  • 

10.  a  N5  =  150 

11.  READ  VALOE  N7  PROMPT'ODTPOT  LINES  (K)  ?  • 

12.  ♦♦)  //R3893XS2  JOB  ( 1 1 35200G,  2G05,  XN  1,  XN7,  %N  5 )  ,  *  JOHN  J  CormACK* 

13.  IF  (N  1  LE  600)  **1  //  H  SGCLASS=Q,  NOTIF  I=R3893 

14.  IF  ( N 1  GT  600)  *♦  1  //  NOTIFY=R3893 

15.  a  N2  =  0 

16.  READ  STRING  S4  PROMPT'PDS  OR  SEQ?  • 

17.  IF  ( *  XS4  *  EQ  ’  PDS  * )  READ  STBING  S3  PROMPT*  DATA  INPUT  MEMBER  NAME 

18.  IF  {'  XS4*  EQ  *  SEQ*  )  READ  STRING  S3  PROMPT»INPUT  DATA  SET  NAME?  ' 

19.  READ  STRING  SI  PRO  HPT  •  GSC24.  F8XX  .  D  AT  A  FILE  NUMBER?  * 

20.  d  SI  =  LPAD  (SI, 2, *0' ) 

21.  I F  (N 2  NE  0)  X  31.000 

22.  READ  STRING  S5  PROMPT'LIST  INPUT  DATA  SETS?  • 

23.  IFCXS5*  NE  •  Y*  )  X  31.000 

24.  **1  //PRINTAXN2  EXEC  PGM =UPRI NT , REGION=60 K 

25.  *«-1  //SYSPRINT  DD  SYSOUT=A 

26.  IF  ('  X  S4  •  EC  *  PDS  '  )  **1  //SYS  UT  1  DD  DSN  =  R3893.  DATA  (  XS3)  ,DISP  =  SHS 

27.  IF  ('  XS4 '  EQ  *SEQ’)  *  +  1  //SYS  OT1  DD  DSN  =  R3893.  XS3,  DISP  =SHR 

28.  **1  //PR INTBXN2  EXEC  PG M= UPRINT, REGION=60K 

29.  *+1  //SYSPBINT  DD  SYSOUT=A 

30.  //SYSUT1  DD  DSN= R3893.  GSC24.  F8XS 1. DATA, DISP=SHR 
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31.  a  12  =  N2  ♦  1 

32.  BEAD  VALUE  HO  PBOHPT'#  OF  PLOTS?  • 

33.  BEAD  VALUE  H6  PBOHPT'#  OF  CTL  FOB  TEKP  DSH?  • 

34.  i  M  =  Hi/60 

35.  a  H8  =  HI  -  60  *  H4 

36.  BEAD  STBIHG  SO  PBOHPT' BIT  BATE  (  INTEGER)  ?  • 

37.  IF  (SIZ  ( SO)  GT  5)  X  36.000 

38.  ♦♦I  //GOXSOXS2  EXEC  PGH=  GASPS  IH,  BEGION=XN  5K,  T  IHE  =  (XN4,XN  8) 

39.  *+1  //STEPLIB  DD  DSN =H38 93 .G A SP. LOAD ,DI SP= SHE 

40.  **1  //FT03FOO 1  DD  * 

41.  BEAD  VALUE  H3  PBOHPT'PLOT  EVERY  XTH  POINT?  • 

42.  BEAD  STBIHG  SO  PBOHPT'POS  OB  NEG  STUFF?  • 

43.  IF  (•  XSO'  EQ  *POS')  a  N1  =  -1 

44.  IF(*  XSO'  EQ  'KEG')  a  N1  =  1 

45.  IF  (  •  %S0  *  EQ  'POS')  X  48.000 

46.  IF  (•  XSO'  BQ  'HEG')  X  48.000 

47.  X  42.000 

48.  *♦  1  XN3,XH1 

49.  **1  /* 

50.  **1  //FT 04  F00 1  DD  DS  H=  B3893.  GSC2  4.  F8XS  1.  DATA,  DISP=SHR 

51.  IF  ('  XS4'  BQ  '  PDS* )  **1  //FT05F001  DD  DSN=B3893.  DATA(XS3)  ,  DIS  P=S  HR, 

52.  IF('XS4'  BQ  'PDS')  *♦  1  //  LABEL=  (,  ,  ,  IN ) 

53.  IF  (•  XS4  •  EQ  '  SEQ  *  )  *+1  //FT05F00  1  DD  DSN=R389  3.XS3,  DISP=SHB 

54.  1  //FT06F00 1  DD  SYSOOT=A 

55.  **1  //FT11F001  DD  DS  N=6GPL0T  1,  DIS  P=  ( N  EH  ,  DELET  E)  ,  SPACE  =  (C  YL  ,%N6)  , 

56.  //  DCB=  (BECFH=VS)  ,  UHIT=SYSDA 

57.  IF  (  NO  EQ  2)  **1  //FT  1  2F001  DD  DSN=SGPL0T2  ,  DIS  P=  (NEH  .DELETE  ),  S  PACE=  (C 

58.  IF  (NO  EQ  2)  **1  //  DCB= (RECPH=VS) , UNIT=SYSDA 

59.  READ  STB ING  SO  PBOHPT* HOBE  BUNS?  • 

60.  IF  (•  XSO'  EQ  'Y')  X  16.000 

61.  **1  // 

62.  LIST  XHO/L  UHN 

63.  BEAD  STRING  SO  PBOHPT* BON  JOB(Y  OR  N) ?  • 

64.  IF  ('XSO'  EQ  'N'J  X  71.000 

65.  IF(  'XSO'  EQ  •  Y •)  X  67.000 

66.  X  63.000 

67.  FREE  DA(LOAD# GSC24.F8XS1  .DATA, DATA) 

68.  SUB  XHO/L  ONN 

69.  JO 

70.  X  4.000 
71  .  SET  ESC 
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STDFF/FILL  PROGRAB  INPOT  DATA 

5,1 

767  800  833  866  433 

145  107  69  32  526 

1000  1000  1000  1000  1000 

1 00001 0000100001 000010000 
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6 ASP  PBOGBid  IBPOT  DATA 

The  iaput  data  consists  of  two  data  aots,  one  is  the  output  data 
froa  tka  ataff/fill  bit  progxaa  aa  shown  in  Appendix  H,  the  other 
is  Hated  below. 

J  J  COBRACK  24  12  27  1979  100000000000011 

2  8  0  0  1  0  25  4  1  200  30  0  10 

1 PHRAX 
2PHHH 
1  PH  OUT 
2PRODB 
3PSLBB 
4  PHBBBR 
SPRTBK 
6PHSBCB 
7FTHK 
arsLErr 
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11 

6  0 

0.0 

1  ODBC 

1 

1 

0.0 

1.0 

2SSYIC 

1 

1 

0.0 

520.0 

3  OOUTPUT 

1 

1 

0.0 

360.0 

4  IS 00 BCE 

2 

2 
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360.0 

5EPRBBB8 

1 

1 
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360.0 

6TPHTBK 
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1 

36.0 

36.0 

0  1.08-7 

1. 

08-8 

1.08-8 

1.08-4 

1  1 

0 

0.0 

1 

1  7813 

5 

9. 589078-3  1 

4  2 

9.6E-3  09.589078-3 

0.0  90.0 

1.8293E-1  .230107 

9.60018-3  09.589078-3 

0.0  90.0 

1.82938*1  .230107 

9.5999B-3  09.589078-3 

0.0  90.0 

1.8293B* 1  .230107 

9.68*3  09.589078-3 

0.0  90.0 

1.82938-1  .230107 

9.58-3 

2887.0  2.0 
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GASP  PBOGBAH  OUTPUT 

This  oatpat  cos sis ts  of  a  tiae  don a in  wavefors  plot  of  several 
paraaeters  is  the  soltipleser*s  saoo thing  loop. 
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